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ABSTRACT 


A system to display images from data generated by the Naval Postgraduate School 
Infrared Search and Target Designation (NPS-IRSTD) system (modified AN/SAR-8 
ADM) was developed using an Intel 80386 CPU based desktop computer as the base 
platform. This computer was enhanced with a Metrabyte PDMA-16 Input/Output board 
to facilitate data transfers and a Data Translation DT2861 framegrabber board for image 
processing. Images are displayed on a 512 by 480 pixel resolution color monitor. 

Digitized NPS-IRSTD data is stored on recording tape, and accessed for imaging 
during playback at reduced speeds. Data is received into the 386 computer through the 
PDMA-16 using Direct Memory Access (DMA) transfers into Random Access Memory 
(RAM). Data is sorted in RAM into the proper format for display, and then moved to 
a frame buffer in the DT2861. FORTRAN programs using Data Translation DT-IRIS 
subroutines were written to manage image display and processing. 

Images produced using the system showed clearly identifiable features from the 
scene observed by the IRSTD. They also provided valuable insight into the IRSTD 
operation. For example, the detector AC-coupling undershoot was clearly evident. 
Additionally, images showed the number of non-operational detectors to be larger than 
previously noted. 

Recommendations are made for system improvements to enhance speed and 


reliability, and for additional study of the imaging properties of the NPS-IRSTD. 
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I. INTRODUCTION 


The AN/SAR-8 is a passive, scanning, infrared detection system originally 
designed to detect the presence of targets of military interest in a maritime 
environment, designate and track detected targets, and display pertinent track 
information on a video monitor. Data was collected and processed to efficiently 
detect the presence of infrared emission sources while strongly discriminating against 
background clutter. Such design parameters led to a system where essentially all 
data underwent some form of processing prior to display. The video displays 
provided information on targets meeting certain criteria while suppressing all other 
scene information. 

The original prototype AN/SAR-8 IRSTD (InfraRed Search and Target 
Designation), the Advanced Development Model (ADM), was obtained by the Naval 
Postgraduate School (NPS) in January, 1985 after testing and evaluation periods 
both afloat and ashore. At NPS the ADM was placed under the control of the 
Naval Academic Center for Infrared Technology (NACIT) to be used for research 
and student instruction. When it arrived at NPS, the ADM was largely in a state of 
disrepair. As part of returning the system to an operational condition some 
modifications have been performed to enhance its research potential. The most 
significant of these modifications was the bypassing of the background suppression 
hardware on one of the two focal plane infrared detector arrays to allow access to 
unprocessed data. As the result of these modifications the system is no longer truly 
an AN/SAR-8, and thus except for historical discussions will be referred to 
throughout this thesis as the NPS-IRSTD. 

The goal of this thesis project was to develop a method whereby the infrared 
scene data generated by the NPS-IRSTD can be displayed on a video monitor, and 
then by applying various image processing techniques, to enhance the images 
produced. The imaging system is now available for a number of useful applications 


including: 


* rapid visual data searches for locating data on interesting scene compositions 
desirable for testing background suppression algorithms. 


« visual analysis of image after clutter rejection algorithms have been applied. 


* continued research in imaging properties of the system and infrared image 
processing techniques. 


* evaluation of image data collection and processing techniques applied for 
possible wider IRST program applications, specifically examining the viability 
of framegrabber technology for shipboard IRSTD system video. 


* instructional demonstrations on infrared sensing and imaging characteristics. 

In the course of working through this project, a number of decisions were 
made that restricted the focus in an attempt to scale it to a reasonable size given 
the time available, and thereby hopefully provide a sufficient percentage of that time 
to be used to study the physics of the imaging process. The most notable of these 
choices came regarding the rate at which the system should be able to produce an 
image. Because of the high data rate output by the NPS-IRSTD and the potential 
for engineering complications arising from hardware idiosyncracies (known and 
unknown), it was decided not to attempt implementation of an imaging system 
relying on data generated at normal system speeds. Instead the image would be 
produced using data recorded on a tape recorder at normal speeds and then played 
back at a much slower rate. 

The hardware used to receive the data from the tape recorder and produce the 
image was an IBM _ compatible computer based around the Intel 80386 
microprocessor chip. This computer was enhanced with a large 320 MB ESDI hard 
disk, an eight megabyte memory expansion, a high speed digital input/output board, 
a VGA display adapter connected to an 800X600 pixel resolution multisync monitor, 
and a Data Translation high speed image processing board possessing four 
megabytes of on board video image storage memory which outputs to a dedicated 
video display monitor. The Data Translation DT2861 board will be referred to as 
the framegrabber. Figure 1 is a block diagram of the NPS-IRSTD PC based 
imaging system. 

To form an image from the data recorded on tape a direct memory access 


(DMA) transfer is employed to move data into the standard random access memory 
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Figure 1 Block Diagram of the NPS-IRSTD Imaging System 


3 


area of the 386 computer. This data is arranged in a format incompatible with the 
display characteristics of the framegrabber, so it is unscrambled into a correctly 
formatted block and then moved into framegrabber frame store memory. Running 
Fortran programs based around Data Translations DT-IRIS subroutines allows image 
data written to frame storage memory to be processed and displayed on the 
dedicated video monitor. 

The image processing and video display functions have been divided into two 
programs, BEGIN and IMAGE, each of which fills a distinct role. BEGIN is 
designed for use prior to viewing images to build tables which define display color 
characteristics. IMAGE is used to control processing and display for specific 
images. 

The remainder of this report is organized as follows: 


¢ Chapter II provides additional information on the history of the IRST program 
and the AN/SAR-8 ADM as well as a description of the system hardware as it 
exists at this time. 


¢ Chapter III comprises the body of the report. The details of the imaging 
process are spelled out, including information on engineering approaches that 
failed to produce acceptable results. This is followed by a description of the 
image processing techniques employed while attempting to enhance the image. 


* Chapter IV presents the results achieved with the imaging system and 
discusses the pertinent features of these results. It also discusses the 
performance results of the imaging system, noting items which require 
attention in order to increase its performance and reliability. 


* Chapter V draws conclusions based on the findings, and then makes 
recommendations regarding areas for future investigation which became evident 
during the course of this project. 


* Appendix A discusses the operation of infrared detectors, with special 
emphasis on the semi-conductor photovoltaic type used in the NPS-IRSTD. 


* Appendix B discusses the imaging properties of the NPS-IRSTD using the 
concepts of linear systems filter theory. 


¢ Appendix C is a collection of the computer program listings used in 
generating and enhancing the video images. 


¢ Appendix D provides detector calibration data. 


II. THE NPS-IRSTD 


A. PROGRAM HISTORY AND BACKGROUND 

The AN/SAR-8 IRSTD system, including the extensive research and 
development program conducted to bring it into existence and support its continued 
evolution, is the result of the Navy’s recognition of the desirability of a modem, 
effective, passive detection system on board its ships. A passive system can be 
developed by exploiting the fact that all objects at temperatures above absolute zero 
(0 K) emit electromagnetic radiation characteristic of the exposed surface 
temperature and emissive properties. Such a passive infrared detection system 
employed on a ship would complement the ship’s other detectors and consequently 
increase the overall flexibility, reliability and effectiveness of the composite 
detection package. It could be operated independently during a period when 
emission control measures (EMCON) are in effect, or it could be operated 
concurrently to provide an additional source of information to verify information 
received from other sources, or for cuing of other detection assets. 

Electromagnetic radiation emitted from an object by virtue of its surface 
temperature and emissivity is a phenomenon that is not dependent on an 
illuminating light source. For targets of interest such as aircraft, missiles, surface 
ships, etc., temperatures range from ambient atmospheric for some external surfaces 
to 1000 °C for the jet exhaust. All objects in this temperature range radiate 
strongly in the infrared region of the electromagnetic spectrum. 

Motivated by the desirability of such a system, the AN/SAR-8 IRSTD was 
developed to detect, evaluate and track sources of infrared radiation. The following 
is a summary of the major milestones, taken from the thesis by Ayers [Ref. 1], in 
the history of the AN/SAR-8 ADM which is now located at NPS: 


* In 1960, the U.S. Navy sponsored the development and testing of the first 
Infrared Electronic Warfare System (IREWS). 


IRST (first test version on the IRSTD) prototype testing took place by SPAR 
in 1969. 


In 1974, the Canadian Navy sponsored the development and testing of the 
Shipboard Passive Surveillance and Detection System (SPSDS). 


The U. S. Navy posted the operational requirement for a passive IRST system 
in 1975. 


A U.S.-Canadian Memorandum of Understanding and a Joint Project 
agreement was signed in 1976 which resulted in the actual contracting and 
construction of the Advanced Development Model(ADM) of the IRST system. 


May of 1978 the IRST system was tested in Halifax, Nova Scotia. 


In June of 1978 the IRST went aboard HMCS ALGONQUIN for at-sea shake 
down and further testing. 


In April of 1979 the ADM was moved to Key West, Florida for further 
testing and program validation. 


Upon completion of the Key West testing, OPTEVFOR recommended full 
scale engineering development of the IRST/IRSTD. The system was then 
transferred to the USS Kinkaid, in September of 1979, for continued testing 
and developmental considerations. 


After the trial the ADM was sent to the Naval Surface Weapons Center 
(NSWC) at Dahlgren, Virginia. 


Full scale engineering development was approved in June of 1981. 


In 1983 a subsequent Joint Agreement was signed by the U.S. and Canada for 
engineering development of the IRSTD system. 


The prime contract was signed with SPAR and General Electric in August, 
1984. 


In January of 1985, the ADM was transferred to NACIT at NPS Monterey, 
California to be used in research. 


In 1986 the cryogenic cooling engine was replaced by a system which uses 
liquid nitrogen. One germanium window was removed from the optical path 
to accomplish the alteration. 


Thesis by Parker December, 1986 recommends realignment of optical system 
to compensate for window removal. [Ref. 2] 


Subsequent to completion of recommended optical realignment, a thesis by 
Ayers September, 1987 studies calibration of detectors. [Ref. 1] 


Bypassing Background Normalizer Unit (BNU) and initial installation of signal 
transmission lines to Spanagel 210 was completed in 1988. 


* January, 1989 six detector data recorded. Problems with data quality 
discovered. 


* August, 1989 data quality problems corrected, full multiplexed data recorded. 

As can be noted from this list, there have been some changes made to the 
system since at its arrival at NPS. The reasons for change were twofold. First, 
when the AN/SAR-8 ADM arrived after its numerous trials, transfers, and storage 
periods it was found to be in an obvious state of disrepair. Second, its usefulness 
for research in its original configuration would have been limited by the background 
normalizing hardware. 

The first major item modified was the detector cooling system. The original 
cooling system was driven by a Cryogenic Technologies cryogenic engine operating 
on a Gifford-McMahon cycle. This engine was complicated, expensive to repair, 
and had a limited operational lifetime. Additionally, a leak was discovered at the 
top of the inner germanium window which was allowing water to condense on the 
detectors. For these reasons it was decided to replace the original system by using 
a pool of liquid nitrogen as the coolant. [Ref. 3] 

In evaluating the optical system after the cooling system was modified Parker 
analyzed the effects of the removal of the inner germanium window, which changed 
the optical path by removing approximately 4 mm of refractive material near the 
focal plane of the system. Using the Super-Oslo optical design software developed 
by Sinclair Optics, INC. he found that a focal shift of the dewar sub-assembly of 
approximately 0.133 + 0.030" toward the reflecting mirror would restore the optical 
performance by reducing the spot size of a point source to less than the 
characteristic size of the individual detectors (0.02"). The corrected configuration 
spot size would be twice the original spot size, but well within the detector size. 
The MTF of the optical system would also be slightly affected; a reduction in 
modulation of approximately 6% from the original configuration can be expected at 
a spatial frequency of 250 lines/inch. [Ref. 2] 

Following the corrections made to the system verifying Parker’s analytical 


work, Ayers and Crittenden [Ref. 1] performed calibrations of detector responses. 


This work was used in making first order corrections to the data used for imaging. 
A summary of the detector calibration results is included as Appendix D. 

The task of restoring the NPS-IRSTD has required a great deal of effort by 
many people; students, staff and faculty. Special recognition is due Mr. W. J. 
Lentz for his work on the system’s electronics. The major modifications were 
complete and the first sets of data were taken in January-February, 1989. While 
taking these initial data sets the detector multiplexer clock was disconnected so that 
only six detectors on each array were sampled. This data showed intermittent, 
unpredictable instabilities when retrieved from the tape recorder for analysis. 
Subsequent data sets taken in March through July, 1989 showed similar 
characteristics. Figure 2 is a plot showing the data irregularities. After numerous 
diagnostic checks and fine tuning modifications, which revealed a variety of 
problems ranging from crosstalk in the data lines in several locations to a bad 
channel in the tape recorder, the system was made fully functional in August, 1989. 
The data used in the video imaging processes for this thesis was collected August 
12, 1989. 


B. CURRENT SYSTEM HARDWARE CONFIGURATION 

The following description of the NPS-IRSTD ADM is intended to give the 
reader a complete overview of how the system works as it is presently configured. 
The emphasis is on the role of each major component and how those components 
interact. 

The NPS-IRSTD is presently installed on the seventh and eighth floors of 
Spanagel Hall at NPS. The scanner assembly and buffer power unit are located on 
the eighth floor, the highest location on NPS grounds , with a 360° view of 
Monterey Bay and the surrounding hills. The analog-to-digital converters and 
system control panels are located in Spanagel 703. The signals generated by the 
NPS-IRSTD are transmitted to Spanagel 210 for tape recorder storage and computer 
analysis. 
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Figure 2 Unmultiplexed Data from the NPS-IRSTD Showing the Dropout 
Characteristic 


1. Scanner Assembly 
The scanner assembly contains the hardware necessary to collect infrared 
radiation in a 360° field of view from 1/2° below the horizon to 10° above, and 
convert this radiation to an electrical signal that can be transmitted to other parts of 
the system for analysis. The scanner assembly houses: 
¢ Schmidt F/1 catadioptric telescope. 


* Two vertical 90 detector arrays operating in the 3-5 wm range. These are 
mounted side-by-side with approximately 1/2° horizontal separation and are 
referred to as the lead and lag arrays. 


¢ A detector cooling system that uses a pool of liquid nitrogen as the coolant. 
¢ Pre-amplifier bandpass filters. 


* Background Normalizer Units, one for each array. The one for the operational 
array has been bypassed. 


¢ Detector signal breakout and multiplexing hardware. 


¢ Slip-ring assemblies for passing the data signal from the rotating head to the 
fixed base. 


¢ Position-in-rotation signal generating hardware. 
Figure 3 shows the scanner assembly and the buffer power unit mounted on the 


roof of Spanagel Hall. 


a. The Telescope 

The Schmidt catadioptric F/l telescope receives infrared radiation 
through a 10" corrector plate. The corrector plate is made of germanium which 
transmits well in the infrared and has the necessary material properties required of 
lenses. The light is focussed onto the detector assembly by reflection from a 16" 
diameter spherical mirror at the back of the telescope. The detector assembly 
consists of two focal plane detector arrays on a cooling finger enclosed in an 
insulated cylinder. Light reaches the detectors through a germanium window 
mounted in the cylinder wall. The two detector arrays are also covered by separate 
filters which selectively pass different portions of the 3-5 [im infrared band. Figure 


4 shows the telescope and detector arrangement. 





Figure 3 Photograph of NPS-IRSTD Scanner, and the Buffer and Power Unit 
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Figure 4 NPS-IRSTD Telescope and Detector Assembly 
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b. Detectors and Detector Noise Reduction 

The NPS-IRSTD’s detectors are Indium Antimonide semi-conductor 
photodiodes, and they function as quantum detectors. That is, they work by 
counting individual incident photons which have energies similar to the band gap 
energy of the semiconductor material. Photon energies are determined by the 
relationship E=hc/A. Detector output is a voltage signal related to the number of 
photons received. Of relevant interest to this thesis is the fact that production of a 
television type video image was not a principal concem in development of the 
ADM. Consequently detector size, type, scan pattern and other imaging properties 
were designed to specifications more conducive to effective detection, categorization, 
and tracking of targets with certain defined infrared signatures, than to detecting 
differentiations in the full range of scene elements. Appendix A provides additional 
information on the workings of the indium antimonide detectors. 

Because of the small band gap energy of the semiconductor material 
there is sufficient thermal energy at normal atmospheric temperatures to excite 
electrons within the material. This produces a background noise voltage level 
capable of masking the detectors’ ability to register the arrival of an infrared 
photon. This thermally generated noise is controlled by cooling the detectors to 
near 85 K with liquid nitrogen. 

The ability to detect incident radiation is also affected by "1/f" 
noise, a category including a number of mechanisms for which the fluctuation 
amplitude varies approximately reciprocally with frequency. The frequency here is 
not that of the radiation, but is the frequency of the electronic signal. Each of the 
180 detectors is coupled to an electronic channel. The effect of 1/f noise is limited 
by connecting each channel to a pre-amplifier band pass filter to eliminate the low 


frequency components of the voltage signals coming out of the detectors. 
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c. Detector Cooling System 
The modified detector cooling system uses liquid nitrogen filling the 
detector support stalk from a reservoir to cool the detectors to approximately 85 K. 
This arrangement is shown in Figure 5. Other items of interest in this figure are: 


* the regions of foam insulation which restrict thermal transfer between the stalk 
and surrounding areas. This increases the time it takes to boil off the supply 
of liquid nitrogen, thus increasing operating time between refills, and also 
reduces the degree of cooling experienced by surrounding areas diminishing 
the risk of condensation of moisture on easily damaged surfaces. 


¢ the heating elements surrounding the germanium window. These are designed 
to prevent condensation on the window surface that would block the infrared 
radiation and could possibly damage the window. 


* the vent tubes which serve to direct the evaporating nitrogen into the area 
around the detectors and eventually into the telescope enclosure to assist in 
maintaining a moisture free atmosphere. 


d. Signal Multiplexing Hardware 

On the operational (lead) array the signal coming out of the pre- 
amplifier bandpass filter has been re-routed. Instead of connecting to the 
Background Normalizer Unit, where all data except that which contains patterns 
meeting criteria indicative of the presence of a pre-defined target type is suppressed, 
each detector channel is brought into a breakout box to provide direct access to the 
signal. This breakout box contains 90 connection points, one for each detector in 
the array, with the continuously varying detector output voltage at each point. 

To reduce the number of cables running from the scanner to 
Spanagel 703, and to reduce the number of electronic components required to handle 
the data, the 90 detector output channels are connected to six multiplexers, fifteen 
channels to each multiplexer. Each multiplexer selects one of its 15 input lines for 
output. At any time, the six multiplexer outputs are for detectors physically 


separated by 14 detectors. Detector multiplexing is shown in Figure 6. 
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Figure 5 Simplified Diagram of the NPS-IRSTD’s Modified Detector Cooling 
Apparatus 
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Figure 6 Multiplexing Arrangement for the Lead Detector Array 
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e.  Slip-ring Assembly 

To get the electrical signals from the continuously rotating head to 
the fixed base a set of slip rings is employed. There is one slip ring for each 
multiplexer channel, 12 in all, and several others for synchronization/timing signals. 
Each slip-ring assembly consists of a continuous, circular copper contact strip on the 
bottom side of the rotating head. Attached to the fixed base is a contact brush 
which rests in constant contact with the rotating ring. An area currently being 
studied for possible improvement of the NPS-IRSTD is to replace the slip-rings with 


an air-gap coupled fiber optic system. 


f.  Position-in-Rotation Signals 

The NPS-IRSTD generates the signals necessary for correlating scene 
data with its corresponding azimuthal position. This is done using gear driven 
signal generators which produce a certain number of pulses based on gear ratios 
during each revolution of the scan head. Presently there are two position-in-rotation 
signals available for use. The first is 1/60 (changeable to 1/120) of a revolution 
which divides the scan into 6° elements referred to as sectors. The second is 
1/1000 of a sector. The rotation time of the head is usually set as close to 2 
seconds as possible so the 1/1000 signal equates to approximately 60 KHz and is 
generally referred to by that frequency. A single, fixed position pulse has been 
superimposed on the channel carrying the sector pulses to give an end-of-rotation 
mark. This is not generated by a gear driven wheel, but rather by passing a 
metallic tab affixed to the scanner through a notch attached to the base, thereby 


interrupting a circuit in the notch once each rotation. 


3. Buffer and Power Unit 
The Buffer and Power Unit acts as an intermediary between the scanner 
assembly and the electronic processing hardware in Spanagel 703. Because of the 
large separation, approximately 150’ of cable run, between the scanner and the 
equipment in room 703 it is necessary to boost voltage levels and otherwise modify 
the transmitted signals to match the characteristics required at each end. The Buffer 


and Power Unit can be seen in Figure 3 to the right of the scanner assembly. 
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4.  Spanagel 703 Electronics 

The electronic equipment located in Spanagel 703 includes the electronic 
power supplies, the power control switches and associated status indicators, system 
timing and synchronization signal generators and data signal analog-to-digital 
converters and multiplexers, which can be seen in Figure 7. The original ADM 
data processing and display units which interfaced with the existing equipment for 
system testing and evaluation were not suitable for research applications and were 
not installed. 

The components in room 703 for which a functional understanding is 
essential for processing the data generated by the NPS-IRSTD include the system 
timing and synchronization signal generators and the data signal Analog-to-Digital 
(A/D) converters and multiplexer. The fundamental clock frequency for data 
produced by the NPS-IRSTD is 5.4 MHz; that is 5.4 x10° detector digitizations are 
performed each second. All frequencies driving components involved in data 
digitization are derived from this rate. There are twelve eight bit A/D converters, 
two sets of six, with one set for each detector array. The twelve A/D converters all 
connect to eight multiplexers, one for each data bit. It is the output from these 
multiplexers, along with the synchronization and rotation signals, that is transmitted 


to the tape recorder in Spanagel 210. 


5. Ampex Tape Recorder 

The tape recorder used in this project is an Ampex HBR 3000i which 
can record twelve independent channels of data simultaneously. Tape recorder 
operation may be driven by either an internal or an external clock, and tape speed 
is selectable from seven possible choices ranging from 1 7/8 inches per second (ips) 
to 120 ips. Data is recorded at 120 ips using the 5.4 MHz IRSTD clock, and 
played back for imaging at 1 7/8 ips using the internal tape recorder clock set to 
4.32 MHz. 
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Figure 7 Photograph of Spanagel Hall, Room 703, Electronics Racks 
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6. Data and Timing Signal Interrelationships 

This section steps somewhat beyond the stated framework of an overview 
but this is necessary to develop an adequate understanding of the system. The 
NPS-IRSTD 5.4 MHz digitization rate includes all 180 detectors in each period of a 
30 KHz clock, (this 30 KHz is not in any way related to or synchronized with the 
30 KHz rotation signal). Therefore, on the transmission line to the tape recorder a 
digitized data byte from a given detector will be followed by the next byte from the 
same detector 33.3 pis later, with information from the other 179 detector’s 
transmitted in between the two. However, because of signal multiplexing, the order 
in which the detectors are digitized is not sequential from top to bottom. 

To relate the detector digitization sequence these conventions will be 
followed: The detectors on the lead array are numbered 1-90 with #1 on top and 
#90 on the bottom; the detectors on the lag array are numbered 91-180 with #91 on 
top and #180 on the bottom. Starting at a point in the cycle where the top detector 
on the lead array, #1, is being digitized, the sequence continues as follows: 


¢ 1,16,31,46,61,76,91,106,121,136,151,166,2,17,32,47,62,77,92,107,122,137,152, 
167,3,18,33,48,63,78 93,108,123 ,138,153,168,4,19,34,49,64,79,94, 109,124,139, 
154,169,5,20,35,50,65,80,95,110,125,140,155,170,6,21,36,51,66,81,96,111,126, 
128,143,158,173,9,24,39,54,69,84,99,114,129,144,159,174,10,25,40,55,70,85,100, 
115,130,145,160,175,11,26,41,56,71,86,101,116,131,146,161,176,12,27,42,57,72, 
87,102,117,132,147,162,177,13,28,43,58,73,88,103,118,133,148,163,178,14,29,44, 
59,74,89,104,119,134,149,164,179,15,30,45,60,75,90,105,120,135,150,165,180. 


Notice that blocks of six detectors from the lead array (e.g. 1,16,31,46,61, and 76) 
alternate with blocks of six detectors from the lag array (e.g. 91,106,121,136,151, 
and 176). 

To decipher how this sequence is generated, start by looking at the 15-to- 
1 multiplexed channels in the scanner. By setting the multiplexer input address 
lines for line zero, the data from detectors 1,16,31,46,61,76,91,106,121,136,151,166 
is presented on the multiplexer outputs. These analog data signals are connected to 
the inputs of the A/D converters, which generate digitized 8-bit bytes. Cycling 
through the 12-to-1 multiplexer attached to the A/D converters sequentially 
introduces the digitized data for detectors 1,16,31,46,61,76,91,106,121,136,151,166 


onto the line leading to the tape recorder. The sequence for detectors 
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2,17,32,47,62,77,92,107,122,137,152,167 is generated by setting the input address 
lines for line one. This process is repeated until the input address lines reach line 
fourteen, after which they are reset to zero and then the process starts again. 
Figure 8 shows the digitization multiplexing scheme. 

The next complication arises from the signal timing for the various parts 
of the system. All clock signals originate in room 703; however, they do not all 
propagate the same distance to the various devices which require synchronization. 
For example, the rising edge of the 30 KHz digitization signal is used to set the 
scanner multiplexers output selector bits to zero. But to do this, the signal must 
propagate from room 703 to the scanner and reset the multiplexer. The data signal 
on the multiplexer output line must then propagate back down to the A/D 
converters. Assuming signals propagate at the speed of light in both directions, 
given the approximately 300’ of wiring in a round trip and the 5.4 MHz digitization 
rate, nearly two digitization periods will have passed before the signals reach the 
A/D converters. Added to this are the delays encountered in propagation through 
the electronic circuitry. Since all these delays are subject to fluctuations based on 
temperature and on the condition of the wiring, circuit components and connections, 
it is clear that the precise synchronization of signals is a complex issue. 

Having given a flavor of the intricacies of the timing signals, the issue is 
to find a way of sorting out the relationships between these signals to make them 
useful for data analysis. The first step towards making signals usable involves 
synchronizing them in phase-lock-loops where a feedback process constantly adjusts 
the timing to establish fixed phase relationships. Once the signals have stable phase 
relationships it is possible to determine the actual relationships between different 
signals as they are transmitted to the tape recorder by displaying them on an 


oscilloscope. 
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Figure 8 Analog-To-Digital Converter Multiplexing 
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il. EXPERIMENTAL/ENGINEERING PROCEDURES 


A. GENERAL 

Given the desired goal of producing video images of the scene being observed 
by the NPS-IRSTD, this project naturally divided itself into two distinct parts. The 
first task was to develop a method by which the digitized data could be manipulated 
into a form conducive to display on some type of computer video monitor. Once 
this method was fully developed and verified, the next task was to acquire real 
image data and to perform various image processing techniques to produce the best 


possible image. 


B. IMAGING SYSTEM HARDWARE AND SOFTWARE ENGINEERING 


1. Introduction 

Before commencing with the development of the data transfer and 
manipulation processes, several fundamental choices had to be made. The first 
issue related to the speed at which the system should be designed to operate. It 
would be desirable eventually to have a display capable of producing images at the 
same rate as the NPS-IRSTD scans (referred to locally as ’real time’). Considering 
the large amount of data to be handled and the high 5.4 MHz digitization rate of 
the NPS-IRSTD, the goal of real time operation was beyond the scope of this 
project. Instead of trying to handle data in real time, the imaging process would be 
accomplished using data recorded on a high speed tape recorder and played back at 
much slower rates. 

The next issue was the choice of the computer to serve as the first stage 
receiver of data from the tape recorder. The choice here was between two 
completely different computers, a MASSCOMP (a minicomputer built around a 
Motorola 86000 CPU) and an IBM compatible 80386 based microcomputer running 
at 20 MHz (Everex Step 386/20). While the MASSCOMP possessed a number of 


desirable features, including a proven interface with the tape recorder and a large 
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library of sophisticated data analysis software, its use also carried definite 
drawbacks. Because the MASSCOMP video display is limited in capability when 
compared with the two displays available with the 386, a second transfer from the 
MASSCOMP to the 386 would be necessary to take advantage of the 386’s superior 
video monitors. In addition to the obvious awkwardness and slow operation of such 
an arrangement there was the secondary consideration of needing to learn two 
completely different sets of computer operating systems and software applications 
packages before really beginning to work effectively. Consequently, the decision 
was made to commit to finding a valid method of passing the data directly from the 
tape recorder to the 386. 

The final fundamental task was to select one of the two monitors 
available to the 386 as the principle video display for the NPS-IRSTD data. Here 
the choice was between the standard computer video display, an 800 column by 600 
row resolution multisync VGA monitor, and a specialized video display system add- 
on consisting of a Data translation framegrabber board and dedicated display 
monitor. The framegrabber system is specially designed for high speed image 
processing, and thus was the more desirable choice if it could be effectively 


employed. 


2. The Everex 80386 Computer System and Associated Components 
Before getting to the description of the image display process at its 
current stage of maturity, I believe it necessary and beneficial to include a 
reasonably detailed description of the 386 computer system which includes all the 
relevant add-on components. This description is intended to provide a clear 


understanding of the operation of each major component, and how they all interact. 


a. Base 386 Computer, Memory/Mass Storage, DMA controller 
The Everex Step 386/20 is a fully IBM PC compatible desktop 
computer based on the Intel 80386 Central Processing Unit with a maximum system 
clock speed of 20 MHz. In addition to the digital input/output board and the image 
processing board discussed in the following sections, the computer used for this 


project was equipped with: 
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* One megabyte (MB) of system RAM. RAM, Random Access Memory, is the 
workspace used by the computer to mn application programs. A _ standard 
ASCII character requires one byte of memory, and codes employed by 
programs to execute functions or perform manipulations normally require one 
to several bytes per single operation. One megabyte equals 1,048,576 
bytes=(1024 bytes)x(1024 bytes). The low 640 kilobytes (KB) is used for the 
Disk Operating System (DOS) and programs. DOS acts as a translator 
between the operator and the computer to allow programs to be executed, files 
to be accessed, and printers and other peripherals controlled. The RAM above 
640 KB is used for specialized functions such as video display memory. 
Figure 9 shows the system RAM apportionment along with the expanded 
memory and image processor board memory. Note that addresses are given in 
Hexadecimal format. 


¢ Eight megabytes of expanded RAM managed by the Quarterdeck Expanded 
Memory Manager system. Many programs operated by themselves can’t take 
advantage of expanded memory because they have been designed to be 
compatible with earlier computers. Special handlers, such as Microsoft 
Windows or DESQview must be used to take advantage of the extra space. 
Unfortunately, these programs often remap addresses to assign programs to the 
space available and then keep track of them. This remapping conflicts with 
efforts to access specific memory locations directly. 


¢ 350 MB hard drive. Each image stored in a data file requires 1/4 of a 
megabyte. 


* 40 MB tape backup system. This allows fast, regular backup of both data and 
programs stored in case the hard disk system fails. 


¢ VGA adapter and multisync monitor with 800 column by 600 row maximum 
resolution. This display, along with Ontar Corporation’s PC-IMG software, 
allows limited image analysis and processing, in addition to that done on the 
image processing board. 


¢ Intel 8237 D.M.A. controller. This is a specialized circuit, commonly 
available on many computers, but hardly ever used. The DMA controller can 
take over the memory address lines to perform external read/writes directly to 
or from the lower 640 KB of system ram. It is controlled by writing to 
input/output registers 8H (the command register), BH (the mode register), AH 
and FH (mask registers), 2H or°6H (address registers), 3H or 7H (the byte 
count registers), 81H or 82H or 83H (the page register). The data written to 
these registers determines all the characteristics of the transfer to take place; 
the channel used, input or output, increasing addresses or decreasing, number 
of transfers, etc. To function properly the DMA controller must interface with 
a device which sends and receives timing signals in a certain pattern, called a 
handshaking protocol. The input/output board Metrabyte PDMA-16 serves this 
purpose. The PDMA-16 manual provides a detailed description of the DMA 
controller. [Ref. 4] 
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Figure 9 Random Access Memory Organization for the 386 Computer 
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¢ Various input/output ports and registers to interface with peripheral devices 
such as printers, modems, mice, video displays, and disk drives. Certain port 
addresses are to specific types of devices by industry standards, while others 
are left open for use by add-on boards or devices. When installing boards or 
devices, it is often necessary to set a number of switches and run an 
installation program to tell the computer that the board exists at a certain 
location and tell the board to respond to commands addressed to that location. 
To access the devices attached to a port the computer writes command bytes 
to the register/s associated with that particular port. Table 1 shows the 
assignment of ports in the Everex 386. 


TABLE 1 
I/O PORT ASSIGNMENTS FOR 386 COMPUTER 
Address (hex) Device Address (hex) Device 
000-1FF Internal System SUES) TAT 1ejPacL 
200-20F Game 380-38C SDLC comm. 
ZL ZAL 7) Expansion Unit 380-389 Binary comm. 2 
220-24F Reserved 3A0-3A9 Binary comm. 1 
250-26F DT2861 3B0-3BF Mono dsp/LPT1: 
278-27F Reserved aCo=3Ch Reserved 
2F0-2F7 PDMA-16 3D0-3DF Color Graphics 
2F8-2FF COM2: 3F0-3F7 Floppy Disk 
300-31F Prototype Card 3F8-3FF COM1 : 
320-32F Hard Disk 


b. Digital Input-Output Board 
The digital input-output (I/O) board serves as the interface between 
the timing and data signals coming from the tape recorder and the computer. The 
I/O board used for this project was the PDMA-16 which is capable of two basic 
types of input-output modes, programmed I/O and Direct Memory Access (DMA) 
transfers. The PDMA-16 uses a standard 37-pin connector. 16 Data lines are 
divided into two eight bit input-output ports, A and B, and the remaining lines are 
for control and timing signals. Figure 10 shows the pin assignments for the 
PDMA-16. The functions of the pins are as follows: 
¢ AQOQ-A7; Port A data (input/output) 
¢ BOQ-B7; Port B data (input/output) 


¢ A DIR; Port A direction (output only). 0 = Port A is set for input, 1 = Port 
A is set for output. 
¢ B DIR; Port B direction (output only). O = Port B is set for input, 1 = Port 


B is set for output. 


* XFER.REQ.; A positive edge on this input initiates a DMA transfer as long as 
all other settings necessary for transfer are correct. 
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REAR VIEW OF PDMA-16 1/0 CONNECTOR 
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Figure 10 PDMA-16 Input/Output Board Pin Assignments 
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* XFER. ACK.; On receipt of the XFER. REQ. the XFER. ACK. goes low. 
After the transfer is complete, XFER. ACK. retums high indicating valid data 
is on the port in output mode, or that the port is ready to receive new data in 
input mode. 


¢ TIMER OUT; A positive pulse or square wave output from the timer. 


¢ TIMER GATE; this provides a means of inhibiting external and internal pulses 
from the timer. Pulses inhibited when pin is low. 


* INTERRUPT; A positive or negative edge triggered interrupt input selected by 
the interrupt control register. 


¢ AUX1-AUX3; General purpose control output lines corresponding to spare bits 
in the DMA and interrupt control registers. 


¢ +5V; Logic supply from the computer. While it can potentially be used to 
power external peripherals, because it comes directly from the computer power 
supply there is a risk of damaging the computer if used. For this reason the 
lines in the connecting cable for these pins were cut. 


¢ GND; Logic and power return ground. 

Programmed I/O involves reading the data present at the port using 
a software call and writing the data into a variable in the program that is running. 
The data is then available for treatment as would be any other variable used by the 
program. It could be stored, printed, or used in mathematical operations. 
Programmed I/O is generally much slower than DMA transfers. 

DMA writes the data into the standard RAM of the computer 
according to an absolute address. The RAM is divided into 64 KB blocks, called 
pages in the language of the DMA Controller. Starting at OH, page boundaries are 
set every 10000H up to the AOQOOQOH address. The DMA controller can address 
only within a page, thus the maximum number of transfers that can be done at one 
time is 64 KB, and this only when the starting address for the transfer is the lower 
boundary of the page. The timing of a DMA transfer can be controlled by a 
counter on board the PDMA-16, or by timing signals generated external to the 
computer system. 

Setup and control for the various functions performed by the 
PDMA-16 is done through eight addressable registers accessed starting at a base 
address of 2FOH. The addresses are mapped as shown in Table 2. The first four 
addresses correspond to the digital I/O and control, and the last four addresses 
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correspond to the 8254 timer. The 8254 timer is available on the PDMA-16 to 
control the speed of input-output operations, time delays, and interrupts as necessary. 


It was not used during this project. 


TABLE 2 
PDMA-16 REGISTER ADDRESSES 


Address Function Type 


2F0H A port Read/Write 
2F1H B port Read/Write 
2F2H DMA Control Read/Write 
2F3H Interrupt Control Read/Write 
2F4H Counter 0 Read/Write 
2F5H Counter 1 Read/Write 
2F6H Counter 2 Read/Write 
2F7H Control Write 
Status Read 


Because of the variety of operations performed by the PDMA-16, a 
complete overview on its use is not appropriate here. The PDMA-16 comes with a 
manual [Ref. 4] which should be referred to for specifics beyond that provided here. 
Any functions used to transfer data for imaging are thoroughly explained in section 
I1I.B.4, or in the program documentation in Appendix C. Additionally, Metrabyte 
included a set of basic language routines with the PDMA-16 which can be used to 


learn the capabilities of the board. 


c. Image Processing Board and Display 
The image processing board used in this project was the Data 
Translation DT2861-60Hz Arithmetic Framegrabber. The DT2861 has _ four 
megabytes of RAM dedicated to video image storage, a variety of input-output 
methods to handle different video formats, and specialized hardware for high speed 
application of numerous image processing routines. The board can be controlled by 
using one of these methods: 


¢ Data Translations’ IRIS-tutor, an interactive, command line driven software 
package. [Ref. 5] 


¢ Data Translations’ DT-IRIS FORTRAN callable routines. [Ref. 6] 
¢ Direct program writes to various addressable process control registers. [Ref. 7] 
The memory is divided into sixteen frame buffers, 256 KB each, for 


images of 512 rows by 512 columns. However the last 32 rows are not displayed, 


30 


giving a resolution of 512 X 480. Data is stored in memory sequentially by row. 
The base address of each frame buffer stores the data for pixel 0,0 (the upper left 
corner) and the next 511 bytes are row 1 data in order. The address base+(512*n) 
is for the left most pixel in the n-1 row. The data in any one of the frame buffers 
may be output to the video display at a given time. Switching the display between 
one frame buffer and another is done in 1/30 of a second. [Ref. 7] 

The DT2861-60Hz is designed to accept standard video input signals 
including RS-170, RS-330, and NTSC as well as certain non-standard inputs in a 
slow scan mode. These are all analog signals, which are digitized internally for 
frame storage in memory. With an undetermined amount of effort on an interface 
design, it can also be made to accept previously digitized signals. For several 
reasons none of these input methods proved satisfactory for receiving data directly 
from the tape recorder. A more precise explanation of why this was so, and a 
thorough description of how data does get to frame storage is provided in section 
10.B.3. (Ref. 7] 

The DT2861 displays a color or gray level in each pixel on the 
monitor based on the value of the corresponding data byte in memory. The 8-bit 
byte means data values range from 0 to 255. Each value 0 to 255 is assigned a 
specific color or gray level determined by eight bit values set in three separate look- 
up tables; red, green, and blue. For each pixel there are separate red, green and 
blue phosphor spots. The color on the display as seen by the eye is the result of 
mixing the three brightness levels. Table 3 shows examples of output colors. Note 
that gray levels are obtained when all three look-up table values are equal. [Ref. 6] 

In addition to those already described, the DT2861 framegrabber 
board and Data Translation software has the following capabilities: 


*« Create windows, areas within a given frame that can be operated on 
independently. 


* Perform a variety of arithmetic functions including setting all values in an 
active region to a constant, adding/subtracting one frame buffer to/from 
another, and multiplying or dividing all values in an active region by a 
constant. 
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* Perform logical operations including AND, OR and XOR of all values in a 
given frame with a constant and AND, OR, and XOR of one frame with 
another. 


¢ Apply resultant look-up-tables to the results of frame addition to map each 
data value, which range from 0-511, to any new value 0-255 before writing 
the result in the frame buffer. 


¢ Calculate the numerical sum of all the individual data bytes in a frame or 
create an array for a histogram of the distribution of data values. 


* Perform convolutions on the images, or apply various types of filters to 
enhance edges, or other characteristics. 


* Zoom, Pan and Scroll to inspect details of an image. 


* Create graphic overlays to label, annotate and dress up the image displayed on 
the monitor. 
TABLE 3 
SAMPLE OUTPUT LOOKUP TABLE ASSIGNMENTS 


For selected values of an all color display: 


Data Value RED GREEN BLUE Output Color 
50 0 0 255 Blue 
100 0 200 200 Aqua 
150 0 ZN) 0 Green 
200 200 200 0 Yellow 
250 Zo) 0 0 Red 
For selected values of gray scale display: 
Data Value RED GREEN BLUE Output Color 
50 0 0 0 Black 
100 64 64 64 Dark Gray 
150 128 Lis} 128 Medium Gray 
200 196 196 196 Light Gray 
250 255 255 255 White 


3. Imaging System Development 

At the outset of the project it appeared as though there were a relatively 
large number of ways to attack the basic problem of displaying an image from the 
data generated by the NPS-IRSTD. However, as is evident from the descriptions of 
the various components given in the previous sections, limitations in component 
capabilities significantly restrict the number of acceptable solutions. After sorting 
through the options, determining hardware-imposed restrictions and limitations, and 
then looking for a path through or around them, it became apparent that there was 
only one truly acceptable solution without resorting to design and construction of 


additional external hardware interfacing. 


32 


This section provides the information necessary to keep anyone working 
to improve the system from wasting time on processes already found unsuitable. 
The flow chart in Figure 11 provides a way of summarizing the possible data 
transfer methods. The following paragraphs discuss the issues relevant to the 
options shown in Figure 11 for transferring data from the tape recorder to the frame 
buffers for display. 

There appears to be only one method with even the potential for real 
time image acquisition. The DT2861 can transfer data at a 10 MHz rate through an 
external port designed to interface with the DT2858 Auxiliary Frame Processor. It 
may be possible to use this port to transfer data directly from the tape recorder to 
the DT2861 frame store memory. This is probably not a trivial task. 
Experimenting directly with the internal interface of the DT2861 would cause 
notable anxiety for a novice because serious damage could be done to such an 
expensive board in the event of error. 

There are a couple of other points worth making regarding the DT2861 
and real time display. When reading the DT2861 documentation, the 12 MHz 
maximum transfer for the slow scan rate appears appealing. But slow scan is only 
for analog signals. To use slow scan, the data would have to be reconverted to 
analog signals identical to those digitized by the A/D converters in room 703. Even 
if these problems are solved, the data must still be sorted into the proper format 
which is an inherently slow operation because it requires a minimum of 45 KB of 
data that must be individually moved. 

The factors involving the choice of the Everex 386 computer for 
receiving the data from the tape recorder were discussed previously in section 
I0.B.1. While choosing the specific device in the 386 which would interface with 
the tape recorder initially held out the form of offering alternatives, such was not 
the case. The drawbacks weighing against the DT2861 just discussed still apply 
here, and thus need not be repeated. 

Because of the large disparity in projected maximum transfer speeds 
between DMA and programmed I/O (120,000 bytes/second for DMA vs. 5,000 to 
10,000 bytes/second for programmed J/O on a 4.77 MHz PC), a great deal of effort 
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Figure 11 Flow Diagram of Possible Data Transfer Methods 
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and time was spent trying to understand and implement DMA transfers. Speeds are 
quoted for the 4.77 MHz PC because the PDMA-16 is not engineered to take 
advantage of the AT bus [Ref. 4]. Before the operation of the PDMA-16 was well 
understood, it appeared the only way to implement a DMA transfer started by 
rotation timing signals was to design a simple Transistor-to-Transistor Logic (TTL) 
circuit to count the sectors. The TTL circuits proved beneficial in investigating 
PDMA-16 operations using artificial signals generated on the breadboard. In this 
manner it was found that DMA transfer rates of 100 KHz were easily handled. 
Additionally, Port A and Port B were found to be capable of independent operation, 
allowing Port B to be used instead of the TTL circuits for sampling rotation and 
timing signals. 

Progression from this point did not go smoothly. Two independent 
problems appeared and persisted in disrupting attempts to read data from the tape 
recorder. One was an apparent loss of data. Transfers set up for 64 KB resulted in 
only 12 KB of data appearing in RAM. Suspicions, which finally proved to be 
true, pointed to the DESQview multitasking environment’s remapping of system 
memory as the cause of data loss. Operating outside Desqview solved this problem 
but brought on the appearance of the second problem: regular system crashes. 

After innumerable crashes and reboot sequences, it was discovered that if 
the data levels on the DMA port are not stable, a parity error is generated in RAM 
at the location to which the particular byte is assigned. If that space in memory is 
not accessed, the problem can go unnoticed. It is only when sequential passes 
through RAM are made for data sorting that the system will invariably crash. An 
error check property of Desqview apparently had the ability to prevent the crashes, 
but the unlocatable data could not be used. 

This seemingly left programmed I/O as the only remaining option. 
Experimentation with this method quickly showed that programmed I/O can’t be 
done at high enough speeds to work with even the lowest tape recorder playback 
speed. Using PC-TOOLS’ read-from-data-port command IPORTB by itself in a one 
step loop while the tape recorder played back at its slowest speed resulted in nearly, 


but not exactly, matching data rates. Looking at the data received one could trace 
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out the cyclic pattern of the twelve A/D converters and while at times there were 
twelve bytes before the pattem repeated, at other times there were eleven or 
thirteen. These nearly matched rates occurring with the computer reading as fast as 
possible leave no time to search the timing signals and use them to queue the 
execution of IPORTB calls. 

At this point necessity forced a second, even more thorough look at the 
DMA process. Following many frustrating days of investigation, it was discovered 
that cutting the tape speed to 1/80 of real time and reversing the phase of the 
digitization clock produced successful transfers. 

The initial work on sorting the data into proper format consisted of 
implementing the BASIC functions of PEEK and POKE. PEEK was used to read 
the DATA from RAM and then POKE was used to write it to a file according to 
DT2861 data file format specifications. This proved far too slow because BASIC 
(even compiled BASIC) is inherently slow and because writing to a disk file is one 
of the slowest operations possible. It required two minutes to unscramble and write 
a file representing an image of 512 columns by 90 rows. 

The only branch on Figure 11 not tested is that which represents sorting 
data written by a DMA transfer into an array. Initially there was no known way to 
directly access RAM from inside a FORTRAN program as is required to use DT- 
IRIS to move the array into a frame buffer. Then, when PC-TOOLS routines were 
used to access memory directly using FORTRAN, the arrays could not be used 
because the version of PC-TOOLS available could only be compiled with Ryan- 
McFarlan FORTRAN while DT-IRIS required Microsoft FORTRAN for compilation. 
There are other factors which appear to make sorting to an array undesirable in 
current application, including maximum array size allowed by DT-IRIS and doubling 
the memory required by putting a byte in an Integer*2 space. However, for certain 
types of data processing it may be beneficial to switch to a system of sorting into 
arrays immediately after the DMA and then processing the arrays prior to moving 
the data to the frame buffer. 

Having described many things that didn’t work, it now remains to cover 


those which did. This is covered in the next section. Before beginning that it is 
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worth restating the point that while it appears data transfer rates may be increased 
somewhat by correcting certain recognized problems, it will probably require 
extemally constructed hardware interfaces to achieve significant increases in 
performance. One possibility already discussed was the external port on the 
DT2861. Another would be to build a buffer memory bank in hardware in which 
the data is received in a scrambled format using one clock, and read out to the 


computer in the proper format using a separate clock. 
4. Current System Operation 


a. Data Transfer Timing 

The process of timing the start of a set of DMA transfers to select a 
particular region of the 360° scan proved to be more difficult than anticipated. The 
original concept was to build external circuitry utilizing the End-of-Rotation and 
sector pulse data line to keep track of the position in rotation. Further investigation 
of the capabilities of the PDMA-16 board indicated it was feasible and potentially 
more convenient to track rotation by reading the state of the EOR and sector pulse 
line at the I/O port through software programming. 

Counting the sector pulses at the I/O port proved to be unreliable. 
Noise on this data line provides false "sector pulses" which throw off the counting 
process. Checking this line with an oscilloscope seems to indicate the noise may 
not be just crosstalk from the data lines. For this reason, a method was developed 
to keep track of the position in rotation by timing the process using the internal 
clock in the 386 computer. Unfortunately, the rotation speed of the scanner 
assembly fluctuates enough to make this method ineffective in repeatedly pinpointing 
particular positions. Such a capability is essential in accomplishing frame to frame 


subtractions to look at time varying changes in a scene. 


b. Data Transfer Process 
Image data is moved from the tape recorder to the image processing 
board frame buffer in a three step process: 
¢« DMA transfers into standard computer RAM. 


« A byte sorting move to a separate location in standard RAM. 


ay 


* Large blocks of sorted data moved into a frame buffer. 
This process is shown in Figure 12. Once in a frame buffer, the data will remain 
there until it is overwritten, cleared by a software operation, or the system power is 
turned off. Frame display is controlled by software and will be discussed shortly. 

The DMA transfer results in data from all 180 detectors being 
written into RAM in a column oriented format, where 180 contiguous , RAM 
addresses contain data from single vertical scans down the two arrays in a pattern 
described in section I.B.2. The DMA one time transfer limit of 64 KB contains 
space for only 364 successive scans of the operational array. Three DMA transfers, 
(starting at 70000H, 80000H, and 90000H), are performed to provide sufficient data 
to construct two 512 column wide images, each representing 3° of horizontal scan. 
Since the DMA transfer takes several software operations to set up, one vertical 
scan of the detectors must be skipped moving between pages of RAM. 

The next step, sorting the data into the correct format is done 
completely within standard RAM using BLKMOV, a subroutine from ONTAR 
Corporation’s PC-TOOLS library. It would be desirable to be able to treat all RAM 
equally and sort directly to the frame buffer, but most software is limited to 
addressing the first megabyte of RAM in one way, and any extra memory another. 
In PC-TOOLS, the routine BLKMOV will move bytes (up to 64 KB) between any 
two locations in the first megabyte of RAM. This restriction to the first megabyte, 
not noted in the documentation, was verified by a phone call to ONTAR. [Ref. 8] 

A second routine, BLKMAT, allows addressing above the first 
megabyte, but the smallest addressable unit is the 16-bit word. Since any two 
consecutive bytes written by the DMA transfer are moved to widely separated 
addresses during the sort, this routine is useless for unscrambling the data. It is 
also much slower than BLKMOV, requiring approximately 1 ms to complete 
compared with 1 ps. 

BLKMAT is used to move the correctly formatted image into a 
frame buffer. BLKMAT is limited to moving 32 Kilowords (64 KB) per transfer 
making it necessary to do four calls to the BLKMAT routine to transfer the full 


image. With sufficient data available for two 3° scans, two separate sets of sort 
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Figure 12 DMA Transfer and Sorting Sequence 
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and block transfers are necessary. The first 3° is placed in frame buffer 0 (base 
address = AQOOQOOH), the second in frame buffer | (base address = A40000H). 

To utilize as much of the full display as practical and make each 
individual detector more apparent, each detector is represented by five consecutive 
rows of pixels on the display. Thus a full image requires 225 KB of RAM. The 
225 KB immediately below the DMA region, starting at 37COOH, is used to build 
the correctly formatted image. Table 4 is a printout of selected portions of RAM 


showing how the data looks before and after sorting. 


C. IMAGE DISPLAY AND PROCESSING 


1. Objectives 
The previous sections have illustrated the mechanisms involved in 

transforming the scene information incident on the detectors into digitized data 
stored in frame buffers in the correct format for display and processing. Once this 
task is complete, there exist several possible approaches to processing the data, 
dependent on the perspective of a particular researcher or system user. Working 
from the premise that an IRSTD is primarily a target detection and tracking device, 
many common data processing methods, such as the AN/SAR-8 ADM’s BNU, are 
designed to identify regions or values within the data which exceed some 
predetermined background threshold. In the academic applications envisioned for 
the NPS-IRSTD, the evaluation of background clutter and noise rejection algorithms 
is only one of the possible avenues of research. Other areas of interest include: 

¢ Evaluation of the infrared background scene characteristics 

¢ Testing of advanced hardware modifications 


¢ Use as an instructional platform for study of general IRSTD characteristics, 
such as optical performance, detector performance, and hardware operations 


With these applications in mind, a major objective guiding the 
development of the imaging system was to create a framework which allows study 
of the modifications introduced into the data by the system. Implied here is the 
desirability of eventually removing the system effects to restore the data to the 


closest possible representation of the actual thermal distribution present in the scene. 
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TABLE 4 
SELECTED RAM LOCATIONS AFTER DMA TRANSFERS AND SORTING 
SORTED UNSORTED 


BASE> 40000H 50000H 60000H 70000H 80000H 90000H 


0 126 130 124 Lay 126 134 

al 126 130 27 0 0 0 

2 RAI) 130 127 0 0 0 

3 126 130 WAS) 0 0 0 

4 nee 12s) 121 0 0 0 

ro) 126 130 124 0 0 0 

6 126 130 126 0 0 0 

U 126 29 25 126 126 126 

8 126 128 124 128 126 37 

9 126 128 124 129 130 130 
10 126 U29) LAS 126 EG AS) 
ai LAG 2s) LAS) Lay 125 126 
i ay Las) 124 25 LAD Ws 
its) Lary Ze) 125 0 0 0 
14 126 130 LAY 0 0 0 
AS) ZAG 131 LAS 0 0 0 
16 126 128) 124 0 0 0 
1.7) 126 128 126 0 0 0 
18 126 128 Lay 0 0 0 
Le 126 Lag) Las 126 126 126 
20 126 130 124 126 126 126 
aul 126 123) 125 130 130 130 
22 WAG 130 127 126 L2G) 126 
as 126 130 126 126 126 LAG 
24 126 alShil 124 AALS 114 124 
25 126 130 124 0 0 0 
26 126 129) 126 0 0 0 
Bil 126 129) 126 0 0 0 
28 126 129) 124 0 0 0 
Bo a7 126) 124 0 0 0 
30 126 LAs) 126 0 0 0 
Sil 126 130 Lay LAT 126 126 
G2 126 Lae) LZ 126 126 126 
Shs) 126 128 12S) 130 130 sal 
34 126 128) LAY 126 128 LAT 
38 126 130 128 126 12S} 126 
36 126 ial 126 0 0 0 
37 Ley 130 124 0 0 
38 Lay eZ 9) L207 0 0 0 
Sie 126 Lae 128 0 0 0 
40 126 LAs, 2a 0 0 0 
41 126 130 126 0 0 0 
42 e297, Las) 2a 0 0 0 
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Essential to success in achieving this goal is the effective presentation on the 
display of the features determined to be of particular interest. 

Satisfactory accomplishment of these tasks will allow pursuit of other 
aspects of the analysis with greater confidence in the nature and reliability of the 
data. Removing the system’s influence would render the data, and any results 
obtained from its analysis, more valuable to researchers working on other studies 


related to thermal imaging or detection processes. 


2. Techniques and Applications 

Once the task of filling the frame buffers with data is complete, display 
of the data stored in a frame buffer on the video monitor is controlled by values 
written to various read/write registers on the DT2861 which are addressed through 
I/O locations. However, the register level operations are transparent to the system 
user when DT-IRIS_ subroutines are used. These subroutines use standard 
FORTRAN parameter format to pass arguments between the main program and the 
subroutine. DT-IRIS documentation for each subroutine lists the required and 
optional parameters with the appropriate data type and value range to make 
application in programming generally straightforward. 

Two programs have been written, BEGIN and IMAGE, which utilize the 
DT-IRIS subroutines to provide control over a variety of image processing functions. 
The commented source code for these programs is provided as Appendix C. The 
remainder of this section is intended to amplify some of the more significant 
features designed into these programs. 

BEGIN was principally designed for use in setting display color 
characteristics controlled through output lookup tables. The program begins by 
setting the DT2861 into display mode and selecting frame buffer 0 as the buffer for 
display. Additional operations are selected from a menu which includes: 

¢ Create color scales for frame buffers 0 and 1 
¢ Change output lookup tables 2-7 (Gaussian) 
¢ Change output lookup tables 2-7 (Linear) 

* Clear a designated frame buffer 
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« View a designated frame buffer . 

Since the 90 individual detector channels, expanded to five rows each, 
utilize only 450 of the 480 rows displayed on the screen, the remaining 30 rows are 
used to form an on-screen, full-time color/gray scale reference. These rows are 
filled, in two column pairs, with increasing values from left to right, such that 
columns 0 and 1 contain the value 0 while columns 510 and 511 contain the value 
of 255. The values in a given region are then identified using DT-IRIS graphic 
overlay capability to print the corresponding numeric value onto the screen. The 
intent of this scheme is to provide an instant reference of the value represented by 
any pixel in the display. 

As introduced in section III.B.2.c, the color displayed by any given pixel 
is the result of mixing the relative brightness of the red, green, and blue values 
assigned in the active output lookup table for the data value represented by the 
pixel. The DT2861 has eight output lookup tables numbered from 0 to 7, all of 
which are initialized automatically the first time a DT-IRIS program is run. 
However, any of these buffers may be modified in any manner by the user using 
DT-IRIS subroutines. The two choices on the menu in BEGIN provide two 
different methods for reassigning the characteristics of any of the two output lookup 
tables 2-7. A strong factor in the design of the output lookup table changing 
routines was the desirability of being able to define all the values in the table with 
a minimum number of parameters. Output lookup tables 0 and 1 are purposely 
excluded from changes because of their general applicability as monotonically 
increasing and decreasing gray scales respectively. 

Two terms are important when discussing assignment of values to lookup 
tables: the index and the value. The index refers to the data value the pixel will be 
representing, while the value is the relative brightness to be associated with that 
index. In the "Gaussian" option, values are assigned relative to their relationship to 
a designated mean index. The mean index is assigned a brightness value of 255, 
and values fall off to either side according to the input standard deviation. In the 
"Linear" option, each color scale is computed from two lines connecting three points 


designated by index and value. Figure 13 shows the resulting lookup table 
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assignments for the case of a gaussian with red(mean=255,standard deviation=80), 
green(mean=128, standard deviation=80) and blue(mean=40, standard deviation=60). 
Figure 13 also shows a linear set with red(low index=100, low value=0, middle 
index=255, middle value= 255, high index=255, high value=255), green(low 
index=90, low value=0, middle index=130, middle value= 255, high index=170, high 
value=0), and blue(low index=0, low value=30, middle index=100, middle value= 
255, high index=150, high value=0). 

The majority of the control over image display and processing has been 
built into the IMAGE program. The routines in this program are those which 
directly apply to the display and processing on a particular image. Again, when 
this program is run the user is provided with a menu of options, which in this case 
include: 

¢ Select the frame buffer to display 
¢ Select the desired output lookup table 
* Make a 6° image in buffer 2 
¢ Correct an image for detector responsivity differences 
¢ Fill in the lines on the display corresponding to bad detectors 
* Apply MTF corrections 
¢ Smooth out the vertical 1 detector to 5 row expansion 
¢ Store a frame to a file 
* Retrieve a frame from a file 
* copy one frame from another 
¢ Subtract one frame from another 
The first two choices here simply allow the user to view any of the sixteen frame 


buffers using any of the output lookup tables. 


a. Make a 6° Image in Buffer 2 
The third option, make a 6° image in buffer 2, horizontally 
compresses the data which has been written into buffers 0 and 1, which after a 
normal DMA transfer sequence will consist of two adjacent 3° images, into a single 


image and places the result into buffer 2. This is accomplished using the 





Figure 13 Output Lookup Table Assignments. (Top) Gaussian. 
(Bottom) Linear. 
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characteristic array processing features of the DT2861 by dividing all the values in 
frame 0 by 2 and placing the result in frame 2 and dividing all the values in frame 
1 by 2 and placing the result in frame 3 using a single call to the DT-IRIS routine 
ISDIVC for each frame. Then within each row each adjacent even/odd column 
pair are added together to form a row of 128 data values. The resulting row of 
128 from buffer 2 is written to the left side of buffer 2, while the row of 128 from 
buffer 3 is written to the right side. 


b. Responsivity Correction and Row Filling 

In thermal imaging or thermal detection systems variations in the 
response characteristics of the detectors are expected, including a percentage of 
detectors which fail completely. The fourth and fifth choices on the IMAGE menu 
address these deficiencies for the output of the characteristics of the NPS-IRSTD. 
Using the calibration values determined by Ayers [Ref. 1], all detector values are 
adjusted based on the given deviation of the response values from the mean 
response of the detector array. To prevent dead detectors from biasing this process, 
they are assigned a response value equal to the mean response value. The dead 
detectors may be masked by using the routine called by the fifth menu choice 
which simply fills the corresponding rows with values interpolated between the row 


above and the row below the dead detector. 


c. MTF Corrections 

The MTF routines are designed to operate on the data in such a way 
as to remove particular aspects of the system transfer functions from the image. 
Specific routines have been developed for the x-direction detector scan 
deconvolution, the x-direction optical effects and the y-direction optical effects. No 
y-direction detector array function was included because the first zero of this 
function lies at a spatial frequency approximately equal to the minimum spatial 
frequency in the image. Additionally, most of the other frequencies in the y- 
direction detector function are so close to zero that dividing by them for 


deconvolution would unreasonably amplify small fluctuations in the scene. 
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The corrections are accomplished through application of the 
convolution theorem, which allows the system function to be simply divided out 
from the data in spatial frequency coordinates. The frequency spectrum of the 
image is computed by extracting a row (or column) of data, and applying the FFT 
using FORTRAN subroutines which were written by NPS Professor E. A. Milne. 
According to LLoyd [Ref. 9] and others, the x-direction detector function may be 
computed from the Fourier transform of the function representing the detector shape, 
which in the one-dimensional analysis of a rectangular detector produces a transfer 
function of Sinc(af,), where @ is the width of the detector. The x-direction and y- 
direction optical functions are derived from Gaussian approximations to the MTF 


curves shown in Parker’s thesis [Ref. 2]. Thus, 


2 2 
Mppee 


where f is the spatial frequency, and o’ is the variance. The curve fit was 
established by calculating the variance such that the modulation value computed 
with the resulting Gaussian equation for the highest frequency present in the transfer 
array matched the corresponding modulation value from the curves. 

Here again, the program design is aimed toward investigating the 
effects of changes in simple system parameters without major system changes. To 
look at other system functions, the MTF array may be changed by changing the 
defining parameter in the given equation, changing the form of the equation 
completely, or with a bit more effort by assigning the transfer properties of each 


spatial frequency specifically in a data array. 
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IV. RESULTS 


A. GENERAL 

This project provided a unique overall performance measure in the form of the 
images it produces. Measured by this standard the system provides good results; it 
correctly and efficiently receives and sorts the data and displays images which 
distinctly show recognizable and identifiable features present in the scanned scene. 
In addition, the value of the image processing capabilities which can enhance 
features deemed to be of particular interest or which can be used to study general 
operational characteristics of the NPS-IRSTD may be readily appreciated from the 
results of the preliminary usage accomplished to date. 

Having overviewed the positive accomplishments, it seems best now to turn 
attention towards those aspects which did not fully meet all the standards or 
capabilities envisioned for the system while it was in the design stage. This section 
is intended to provide the basis for making the modifications which will 
undoubtedly improve the systems performance and utility. These topics can be 
divided into the following categories: 

¢ Data transfers 
¢ Image processing program performance 


¢ Output image quality and usefulness 


B. DATA TRANSFERS 

As detailed in section I.B.3, the DMA transfer provided the most direct and 
efficient method of making data available for display without engineering a 
completely new interface device. The DMA process does effectively accomplish the 
role it was envisioned to fulfill, however it does so at a maximum rate that is 
inexplicably slow. The current DMA transfer rate is only 67,500 bytes/sec. 
Attempts to play back the data at twice this rate (135,000 bytes/sec) result in the 


cormption of the computer RAM in the region where the DMA transfers are 
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assigned. The subsequent use of the sorting routine to read data and move it about 
results in the computer processor detecting this corruption and crashing the system. 

One additional idiosyncracy superimposed onto this scenario is the apparent 
discrepancy between the I/O board documentation’s description of the handshaking 
protocol expected by the DMA controller and the arrangement which works in 
practice. The top diagram in Figure 14 shows the handshake as depicted in the 
PDMA-16 manual. From this one would expect to find optimal performance with 
the data changes furthest away from the rising edge of the clock pulse as 
represented by the center diagram of Figure 14. However, trying to operate the 
system in this manner causes the type of crashes previously described to occur. 
Operating the system as depicted by the bottom diagram of Figure 14, where data 
signal changes take place immediately after clock changes, is the only way found to 
make the system work. 

Another characteristic feature of the transfer process is the intermittent 
introduction of an extra byte into the data stream such that one cycle through the 
A/D converters contains thirteen bytes instead of the necessary twelve. This has the 
unfortunate effect of fouling the sorting process for all the subsequent bytes in that 
transfer. Such an event produces a blank band in the image and shifts the 
remaining parts of the image out of position as seen in Figure 15. 

This type of image corruption takes place frequently enough to be annoying, 
but not so often as to render the system unusable. The cause may be as simple as 
noise being introduced on the clock line by signals propagating in the data lines of 
the ribbon cable connected to the PDMA-16. Cutting the length of this ribbon 
cable which connects to the I/O board seemed to reduce the error rate noticeably, 


and changing to a shielded coaxial cable may solve the problem altogether. 


C. IMAGE PROCESSING PROGRAM PERFORMANCE 

Because the majority of the routines in the image processing programs BEGIN 
and IMAGE are direct applications of the DT-IRIS subroutines, the overall 
performance can be summed up by saying that almost everything works as designed. 


49 


inpuc 


| 
A 
Ru a a - 50nS min. 






XFER. ACK. 
Cee?) 60nS typo: a 
output i 





Ty, = D.M.A. transfer time 
1.3 to 5 microsec. for bytes 
& to 8 microseconds for words 


PDMA-6 HANDSHAKING PROTOCOL 


Digitization 
Chek 


2 


ae HAH) EAT See 
EXPECTED PHASE RELATIONSHIP 





Digitization 
ee ee 


| 
SO 
Significant 4A 4) BAT Sec 


Data Bit 
Beemer NSE) een 


FUNCTIONAL PHASE RELATIONSHIP 


Figure 14 PDMA-16 Handshaking Protocol. 
(Top) PDMA-16. (Center) Expected Synchronization. 
(Bottom) Functional Sychronization. 


50 


Bee, 
gg ots 


g 
ee 
< 
= 
2) 
4o} 
Gs} 
[28] 
wo 
= 
© 
a 
OD 
5 
3 
A 
vo 
me 
ov 
OD 
3s 
E 


Figure 15 





However, there are several problems which were not anticipated which diminish the 
overall quality of the final product. 

First and most irritating of the problems resulted from an unnoticed design 
characteristic built into the DT2861. While the DT2861 has 16 frame buffers, its 
address bus is only designed to address any two at a time. Thus when certain types 
of operations are performed in the IMAGE program, such as reading an image from 
a data file, upon completion of this process the address bus remains pointed to this 
particular buffer. If the IMAGE program is terminated at this point and the DMA 
transfer program is run, the transfer will proceed to write the data to the last used 
buffer instead of buffer 0 as designed. When such a problem occurs, the quickest 
and most reliable method of resetting the bus to direct DMA transfers correctly to 
frame buffers 0 and 1 is to use the BEGIN program routine to rewrite the color 
scales on the bottom of these buffers. 

The other problem which relates to the DT-IRIS functions is an apparent flaw 
in the operation of the ISPUTP subroutine, the routine which is used to write data 
values directly to one or more consecutive pixels in a designated frame buffer. 
When ISPUTP is used to write more than one pixel, it seems to write meaningless 
values to the next two consecutive data locations. This only becomes a problem 
when a program, such as the smoothing routine in IMAGE, needs to read and use 
values from locations just written to by ISPUTP. 

The remaining program performance comments relate to the MTF routines in 
the IMAGE program. These programs are not inherent functions of the DT2861, 
DT-IRIS system but instead only use DT-IRIS to access the data present in the 
image. As presently designed the MTF routines are subject to: 


¢ Execution times which are noticeably longer than the other processing 
procedures 


* Effects at the edges of the image due to the continuous data assumption of the 
FFT subroutine 


* Modulation effects which raise the data values in regions containing high 
spatial frequency components above the maximum allowed value of 255 in an 
8-bit byte storage location 
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The speed of execution issue is not unexpected considering the number of 
mathematical calculations and data transfers which must be performed. In fact, 
anyone familiar with performing transform calculations on older, slower computers 
may be initially impressed by how rapidly the entire process can be accomplished. 
X-direction MTFs complete in approximately a minute, but Y-direction MTFs take 
around fifteen minutes. The reason for the difference is twofold. The first reason 
is that there are only 90 rows to be processed in the X-direction case, while there 
are 512 columns in the Y-direction. Secondly, the routine to extract data ISGETP 
can be used to retrieve a full row with a single call, while 90 calls are necessary to 
get the data for a column. One positive aspect worth mentioning at this point is 
that at the current speeds of operation, the user can see the results of the transform 
process as it works its way through the image. 

The edge effects noted in item two of the list pose one of the most interesting 
problems to anyone trying to remove the system’s influence from the data. A 
common approach to dealing with this artificial discontinuity is to apply a data 
conditioning algorithm to either partially or completely remove it. Two examples of 
this type of processing, Hamming and Hanning, are contained in the last two 
subroutines on the IMAGE listing in Appendix C. When considering whether or 
not to implement one of these routines the issue arises on whether or not the cost 
of modifications done on the data by the conditioning algorithm is covered by the 
benefits gained through reducing the artificial high frequency characteristics of the 
edges. 

A preliminary application of the Hamming routine while the data was in its 
original form, unscaled and with no "DC offset" removal, damped out all the image 
data except for the center third of the screen. This effect is demonstrated in Figure 
16 where the original signal is a cosine wave with a period of 51.2 points (10 
cycles/ 512 points), while the lower shows the signal after modification by the 
HAMMING routine. 

The last issue listed, modulation effects at high frequency content edges, is the 
result of the limitations imposed by the 8-bit byte storage format. Often in the 


images, the edge between a "hot" object and the background is quite sharp. This 
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Figure 16 The Effect of Hamming on a Cosine Wave Signal 
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region of rapid transition to values which peak (and periodically clip) at 255 is also 
the type of region which will be most strongly affected by removing the system 
function because the process involves division by the small fractional values 
associated with the higher spatial frequencies. The question then arises whether to 
develop a method to normalize the entire scene based on the maximum value after 
the deconvolution or to simply assign 255 to all data which exceed this value as the 
result of processing. 

The final topic of this chapter, image quality and usefulness, is best addressed 
by showing images and explaining what aspect of system performance each portrays. 
The images in Figures 17 through 35 are photographs of the framegrabber display 
taken with a 35 mm camera with a standard 50 mm lens at about two feet distance, 
using a 1/60 second shutter speed and an f/# 2 setting for ASA 125 black and 
white film. The diagonal streaks in the images are artifacts of the photographic 
process; probably reflections off the display face. 

For the display of all infrared scenes the framegrabber was controlled by an 
output lookup table designed to show more detail in the 100 to 150 value range 
than the standard monotonically increasing scale, unless otherwise noted. The scale 
was created using the linear color table routine in the BEGIN program setting ali 
colors with the set of points (80,0), (180,255), (255,255). 

The following list contains the important data concerning each of the images 
and describes the features worth noting: 


¢ Figure 17 is a 3° wide image showing an antenna which is located near the 
southwest comer of the roof area supporting the scanner assembly. This 
antenna is labelled "A" in the photograph displayed in Figure 18. This 
antenna does not have the double armed feature that it appears to have in this 
image. This object doubling is evident in a number of the other images and 
the exact cause has prompted much speculation. 


* Figure 19 is the adjacent 3° region to the one shown in Figure 17 and was 
created as the second frame in the same DMA transfer set. The prominent 
object in this image is the crane which is labelled "B" in Figure 18. 


* Figure 20 is the result of using the routine to combine the two 3° images into 
one 6° image within the IMAGE program. Note the difference in appearance 
as the horizontal dimension is now less distorted in scale relative to the 
vertical dimension. 
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* Figure 21 is another 6° image, this one for the area adjacent to that shown in 
Figure 20. The prominent object here is an antenna (also seen in Figure 22, 
labelled with an "A") used by the Electrical and Computer Engineering 
department in basic radar laboratory work. The double image is again 
apparent. Also, the broad vertical band structure seen in all the images so far 
is the image of four tall poles attached to a pallet that was sitting in the 
middle of the roof area close enough to the scanner assembly to be badly out 
of focus. These poles can also be seen in Figure 18. 


* Figures 23 through 26 are a set of images similar in format to the set of 
images 17, 19, 21 and 22. The area of view for this set is around the highest 
part of the roof of Hermann Hall above its main entrance. Figure 27 is a 
photograph of this area provided for comparison. Figures 23 and 24 are 3° 
images, while Figures 25 and 26 are 6°. An observation worth noting at this 
point is the doubling of images regarding the antennas on the top of the roof. 
There are only two antennas, not the four apparent in the image. At this 
point it is worth suspending the listing to discuss the image doubling feature 
briefly. 


Careful observation of the doubling feature shows it to be range dependent. 
Of the three objects where image doubling is most evident, the antenna in Figure 17 
is the closest object to the scanner, the antenna in Figure 21 is slightly more 
distant, and the roof of Hermann Hall is significantly farther away than either of 
these. Measuring the distance between the object and its double (being careful to 
remember that Figure 21 is a 6° image) shows that the separation decreases with 
range. One of the leading contenders among possible explanations of this 
phenomenon is the telescope focus, and can be understood by referencing the 
diagrams in Figure 28 which have been taken from [Refs. 1, 2 and 4]. The two 
diagrams at the top of this figure are energy distribution representations of a point 
source imaged onto the focal plane. The size of these spots as calculated by the 
Super-Oslo software is less than the detector size of 0.02" by 0.003". The figure at 
the bottom is representative of the findings when the telescope is out of focus. 
Note the relationship of the spot size to the detector size, and how the detector 
support stalk shadow produces two horizontally separated regions of the spot image. 
It is possible, if the out of focus condition exists, that the image doubling results 
from the splitting of a single point into the two sides of its corresponding spot in 
the image plane. 


Now to return to the description of the images: 
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* Figure 29 is the result of bringing the images seen in Figures 25 and 26 into 


a 12° image. The intent of this image is to give a perspective on how things 
look when the horizontal and vertical dimensions are brought into proper scale. 


Figures 30 and 31 are the same image as seen in Figure 21, displayed using 
different output lookup tables. Figure 30 uses the standard monotonically 
increasing gray scale while Figure 31 uses a monotonically decreasing gray 
scale. The intent here is to attempt to give the reader an impression of the 
significance of the choice of output lookup tables. It is unfortunate that 
presentation of some color output lookup tables is not possible in the thesis 
format, because the effect color adds is quite striking. I believe it is a 
reasonable assessment to say that proper design of output lookup tables is the 
most significant factor in optimizing the display to bring out particular scene 
features or characteristics. 


The remaining figures all attempt to relate the effect of some particular type 


of image processing. A number of the images may appear to be essentially 


identical, because of the minimal effect of the image processing operation through 


which they have passed. Reproduction of the photographs in the printing process 


may even degrade image quality to the point where there is no discernable 


difference between certain images which do show slight differences on the display. 


This points out the fact that certain processing operations may prove to be of little 


or no real value after they have been thoroughly evaluated. Again resuming the 


description of the images: 


Figure 32, the same basic image as Figure 19, has been corrected only for 
variations in response between detectors. Application of this routine actually 
seems to make the presence of the individual detectors slightly more 
noticeable, instead of balancing them as was hoped. This may be due to any 
number of factors including the time elapsed since the calibration data has 
been taken or possibly the fact that installation on the roof and the passing of 
the data all the way through the system creates effects not present when the 
detector outputs were measured directly from a breakout box. 


Figure 33 is the result obtained when the bad detector lines have been filled 
in by linear interpolation. Note that the simplistic gradient calculation 
produces some truly undesirable effects at the edges of tilted or horizontal 
objects, such as has occurred to the crane in this image. 


Figure 34 is the result of performing the deconvolution only for the detector 
shape in the x-direction. The most noticeable effects from this routine were 
apparently the result of the discontinuity at the outside edges of the image, 
which are not real or desired. 
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Figure 35 is the result of only deconvolving the x-direction optical effects. 
Here again this particular process had little impact on image appearance. 


Figure 36 shows the appearance after just the smoothing routine has been 
performed. 


Figure 37 shows the image after all the various processes have been 
performed. 


Figure 38 shows a side by side comparison of images, the left side after being 
processed with ali the possible routines and the right with no processing done, 
put together by the combining routine in IMAGE. 


Figure 39 is also an image on which all the various processes have been 
performed. This is a 12° image of the Monterey hillside which has the 
Sheraton Hotel at the bottom and the Presidio at the top. No positive 
identification of an object has been made, however the prominent vertical line 
in the center may be Franklin street, which can also be seen in the photograph 
in Figure 40. 


Figure 41 is the same image without any processing. There is distinctly a 
difference between the two, but it appears debatable as to which is the better. 
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Figure 17 NPS-IRSTD Image (3° Wide) of an Antenna on Spanagel Hall 
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Figure 18 Photograph the Roof of Spanagel Hall Showing the Antenna 
in Figure 17 and the Crane in Figure 19 
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Figure 19 NPS-IRSTD Image (3° Wide) Showing a Crane on Spanagel Hall 
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Figure 20) : 1 Figures 





Figure 21 NPS-IRSTD Image (6° Wide) Showing the Radar Lab Antenna on 
Spanagel Hall 








Figure 22 Photograph Showing the Radar Lab Antenna 


Figure 23 NPS-IRSTD Image (3° Wide) Part of the Roof of Hermann Hall is 
Visible in the Lower Right Comer 
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IRSTD Image (3° Wide) Showing the Roof of Hermann Hall 





Figure 25° NPS-IRSTD Image (6° Wide) Combining Figures 23 & 24 
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Figure 26 NPS-IRSTD Image (6° Wide) of the Bay Visible to the Right of 
Hermann Hall 








Figure 27 Photograph of the Roof of Hermann Hall 
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Figure 28 Diagrams Representing a Point Source Imaged on the Focal Plane. 
(Upper Left) Orignal AN/SAR-8, Spot Diameter is 0.0009". 
(Upper Right) NPS-IRSTD, Spot Diameter is 0.0016". 
(Bottom) Out Of Focus Condition, Detector Size is 0.02" by 0.003". 
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Figure 31 NPS- 





32 NPS-IRSTD Image (3° Wide) Corrected for Detector Response 
Differences 








Figure 33 NPS-IRSTD Image (3° Wide) with Rows for Bad Detectors Filled 
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Figure 34 NPS-IRSTD Image (3° Wide) Corrected for X-direction Detector 
Shape MTF 
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Figure 35 NPS-IRSTD Image (3° Wide) Corrected for X-direction Optical MTF 
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Figure 37 NPS-IRSTD Image (3° Wide) with all Available Corrections Applied 
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Figure 38 NPS-IRSTD Image Comparing a Corrected to an Uncorrected Image. 
(Left) Corrected. (Right) Uncorrected. 
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Figure 39 NPS-IRSTD Image (12° Wide) of the Monterey Hillside with all 
Corrections Applied 








= WAC a 





4() Photograph of the Monterey Hillside 


Figure 
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Figure 41 NPS-IRSTD Image (12° Wide) of the Monterey Hillside with no 
Corrections Applied 
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A. 


Vv. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


In the process of sorting out what types of concluding comments are 


appropriate for a project such as this, and how they might best be presented, it 


became apparent that there are three distinct topics to be addressed: 


The utility of the imaging system in its present academic environment. 


The potential for generalizing the application of framegrabber technology to 
provide real time imaging capabilities for future IRSTD systems. 


The value of the project to me personally. 


Within the context of the academic applications, the development of the 


capability to display the NPS-IRSTD data in image format using a PC based system 


that is efficient and easy to use adds a valuable asset to the overall NPS-IRSTD 


project. This opinion is supported by consideration of the few uses already 


envisioned: 


It allows the data to be searched rapidly by anyone interested in finding (or 
avoiding) certain features in the system’s total field of view. 


The image processing capabilities support the study of linear filter systems 
theory as it applies to this combined optical and electronic system. This may 
be done in the context of thesis research or as a demonstration/laboratory for 
related courses. 


The results of the studies on system effects could be used to restore the data, 
which has passed through this series of optical and electronic filters, to a state 
which more closely resembles the actual irradiance distribution of the scene 
under investigation. 


The ability to visually present data which has been processed with a 
background clutter algorithm in the form of an image may provide additional 
insight into the algorithm’s overall performance and effect. 


Finally, but certainly not least in importance, is the imaging systems potential 
as a tool to evaluate the performance of the NPS-IRSTD system. Simply 
presenting the data in image form has highlighted the fact that there are 
significantly more bad detector channels present than at the time Ayers was 
working on system calibration. The image display of the scanned scenes also 
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exposed the surprising image doubling effect, which will only be fully 
understood after further investigation. 


On the topic of the potential for employing framegrabber technology to display 
images from future IRSTD systems, it seems appropriate to be somewhat more 
guarded in making an appraisal. This reservation stems from the basic fact that 
many of the problems noticed while using the present system will require solutions 
based on areas of technological expertise beyond the scope of this study. 

Now, having said that, I believe it to be a reasonable assessment that 
framegrabber technology deserves additional study and consideration as a potential 
source of imaging capabilities for operational IRSTD systems. It seems reasonable 
to envision a framegrabber system, connected to an IRSTD through an appropriately 
engineered interface, filling a dual role in a shipboard system. The principal 
function in this environment would be real time display of a selected portion of the 
scan. Such a display could be used by bridge watch personnel as an additional aid 
in nighttime and low visibility navigation, and by CIC personnel in validating the 
position of targets detected by the system. Real time display does not seem 
possible without a significant effort being made to engineer an interface directly 
between the framegrabber and the IRSTD. 

The second role aboard ship, similar to the final item listed for academic uses, 
would be as a tool to evaluate system performance. In addition to highlighting 
obvious items such as bad detectors and unusual imaging behavior, diagnostic 
routines could be developed based on standard infrared system performance 
measures to make checks on_ system _ performance. Properly designed and 
implemented, these diagnostics would fill a role similar to sonar calibrations where 
they provide periodic tests of system performance as well as educate the users on 
the real capabilities (and limitations) of the system. 

Finally, I would like to make a few comments about the value of this project 
to my education, if for no other reason than to remind myself in future years of 
some very important points. At the outset I had the naive impression that it would 
be a straightforward task to bring together the necessary pieces of hardware 


according to the methods prescribed by the available documentation. The system 
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would then, operating as designed, receive and display data only modified by filters 
which match the classical descriptions offered in course lectures and notes; filters 
such as the optical system, detectors and electronics. The task would have been to 
determine experimentally the system MTF. 

It has been an invaluable lesson to see all that can be involved in making a 
complex system such as an IRSTD produce reliable data, and how this can impact 
peripheral projects dependant on the system. In general it now seems much more 
appropriate to approach real world tasks with something less than total optimism 


and also to question, but not dampen, such optimism in others. 


B. RECOMMENDATIONS 
It is believed that framegrabber technology does have a valuable role to play 
in the continuing work of the research with the NPS-IRSTD and the potential for a 
role in operational IRSTD systems. Based on this premise, attention given to the 
following items can only improve system performance and usefulness: 
¢ The DMA transfer process needs to be made faster and more reliable. 


* The source of the image doubling effect should be determined, and if possible 
the problem corrected. 


¢ The problems relating to the application of the FFT routines needs to be 
studied to determine optimal methods for scaling and otherwise handing data. 


¢ The problem with the destination of DMA transfers getting redirected by the 
changing address bus of the DT2861 needs to be addressed. The solution may 
be accomplished best by writing directly to the appropriate DT2861 register. 


* Development of a scheme whereby more than 6° of data can be received with 
any one set of DMA transfers would potentially increase the overall usefulness 
of the imaging system. This could be done by making use of the 8 MB of 
expanded RAM. Six or seven 64 KB transfers could be written to standard 
RAM, and then all data moved to expanded RAM. Then, frame by frame, the 
data could be returned to standard RAM for unscrambling prior to the final 
move to a frame buffer. 


¢ The behavior of the DT-IRIS subroutine ISPUTP needs to be investigated to 
remove its adverse behavior in the smoothing routine. 


* The total system transfer function should be experimentally determined, and as 
an additional step it may prove valuable to make a comparison to theoretical 
predictions. 
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After the system transfer function has been determined, it should be used to 
deconvolve system effects from the data. The significance of the effect of the 
system on the data should then be determined by comparing the raw and 
deconvolved data sets. 


The bad detector channels should be checked for the possibility of restoring at 
least some of them to full operation. Also, the gain on the multiplexer 
outputs should be balanced to remove the horizontal band structure noticeable 
on all images. 
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APPENDIX A 


I. INDIUM ANTIMONIDE INFRARED DETECTORS 
A. DETECTOR PERFORMANCE 

The detectors used by the NPS-IRSTD are made of the semiconductor 
compound Indium Antimonide (InSb), and are connected for photovoltaic operation. 
Each detector is 0.02" high and 0.003" wide and has an instantanious field of view 
of 2 milliradians by 0.3 milliradians. The system contains 180 detectors arranged 
into two arrays, each with 90 detectors in a vertical column. These arrays are 
horizontally separated by approximately 1/2°. The full field of view of 360° by 
10.5° is achieved by scanning the optical sub-assembly in a full circle about the 
horizon. 

Infrared detector performance is generally characterized by certain figures of 
merit which have gained acceptance in the scientific and industrial communities. 
Responsivity, Noise Equivolent Power (NEP), and Detectivity are among the most 
commonly used figures of merit. Figure 42, extracted from The Infrared Newsletter, 
Vol 19, No. 4 [Ref. 10] gives definitions for these terms. Figure 43 is taken from 
The Infrared Handbook [Ref. 11] to show example performance data for InSb 
detectors, while Figure 44 was taken from The Electro-Optics Handbook [Ref. 12] 
to give a comparison between InSb and several other detector types. The values for 
typical InSb photovoltaic detector specifications in Table 5 are extracted from 
Photonics Design and Applications Handbook (Ref. 13]. 

TABLE 5 
TYPICAL InSb PHOTOVOLTAIC DETECTOR SPECIFICATIONS 


FORD AEROSPACE on CORP. 


SPECTRAL Pi/2 RESPONSIVITY OPERATING RISE ACTIV 
rae m) (cmH py ) (A/ W) TEMPERATURE PE asec) AREA {mm 
Se Seater omm) 77K 6x10 to ave 


SANTA BARBERA eggs CTR. 
SPECTRAL RESPONSIVITY OPERATING RISE ACTIVE 


RANGE (_m) (cmH we ) (al) TEMPERATURE slats ia AREA (mm*) _ 
1-5 > Ta 0.25, 0.50 
0, 220 
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Responsivity R(A) 


The responsivity of a detector is defined as the ratio of the output signal in volts or amps 
to the incident optical power in watts. R(A) is measured in volts/ watt or amps/watt. 





= v a 
where v = electrical signal 
®, = incident spectral power at wavelength A, W um"! 
AX = small wavelength interval that includes the wavelength A 
Detectivity D* 


D* is the principle figure of merit for single-element detectors. It is defined as the signal- 
to-noise ratio of the detector when 1 W of radiant power is incident on a detector area of 1 
cm? with a bandwidth of 1 Hz. D* is measured in cm Hz!/? Wl. 


aN, A , Af 
os rege rm 


rms i voltage or current 
= rms noise voltage or current 
rms incident optical power, W 
detector area, cm? 

electrical bandwidth, Hz 
noise equivalent power, W 

= responsivity, VW" or amp W" 


where 


fou fl 


AZP Yr eZw 
5 
i] 


Noise Equivalent Power NEP 


The noise equivalent power is the required power incident on the detector to produce a 
signal output equal to the rms noise output. NEP is measured in watts. 


ims 


Re 


NEP = 


where i, = Fms noise current, amps 
current responsivity, amp/W 


Figure 42 Figures of Merit for Single Element Dectector (Extracted from 
The Infrared Newsletter, Vol. 10, No. 1) 
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Ty 777K 


A, = 10 


Ry = 1 to 100 megohms 


-4 Z 


to 10+ cm 





T< 1 psec a Theoretical 
FOV = 60° > 
Background temperature = 300 K < 
R =10° vw! = 
Pp ~~ 
os Mass ee a0 a 
Pp = 
> 
o 
= 
< 
— "A 
eS 
i.) 
2 
a 
my 
% Signal 
and Noise 
10° noe ile 10° 10° a a ae 
Frequency (Hz) Wavelength (um) 
(a) Frequency response of detector signal (6) Spectral response of detector D* (A, 
and noise. 1000, 1). 


Factor 


* 


Typical Improvement 


Improvement of D 





0 20 40 60 80 100 120 140 
FOV (deg) 


(c) Improvement in D* versus FOV for 
typical detectors. 


Note 


+ (1) Resistance is area and bias voltage dependent. 
(2) D* will increase when the FOV is restricted with cold shields or when the spectral extent 
is restricted with cold spectral filters. 


Figure 43 InSb Detector Performance Data at 77 K 


Handbook) (Extracted from The Infrared 
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Figure 44 Spectral Detectivities for Above Average Crystal Detectors at 
Frequency f,, (Extracted from The Electro-Optics Handbook) 
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B. AC-COUPLING 

Infrared imaging often requires detecting small fluctuations in a large 
background signal. The detectors simply register the arrival of all photons within 
their instantaneous field of view. To extract more information from the data 
generated by the detectors, the average value of the background can be subtracted 
out. A common method of implementing background subtraction in IR systems is 
to couple the detector output to a circuit which blocks the low frequency 
components of the signal. This practice is known as DC-blocking or AC-coupling. 
AC-coupling increases the IR apparent thermal contrast and blocks the 1/f noise, 
improving the signal-to-noise ratio. The NPS-IRSTD detector output is AC-coupled, 
with a low frequency cut-on at approximately 100 Hz, to increase its detection 
capabilities. Lloyd [Ref. 6] gives the circuit shown in Figure 45 as an example of 


a circuit which can be used for this purpose. 


AH — 


INPUT 
Ne = FROM 
DETECTORS 


| ea | 


Figure 45 Simple AC-coupling Circuit 
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AC-coupling does have some drawbacks. One is an inability to calibrate the 
system for determining the actual temperature of objects present in the scene. 
Additional drawbacks include droop and undershoot response characteristics. Droop 


occurs as the detector passes from an extended region of low temperature to an 


Pe 


extended region of high temperature. Undershoot occurs as the detector passes from 
a high temperature region to a low temperature region. Both effects are 
demonstrated by the curves in Figure 46. The top set of curves shows signal droop 
while the center set shows undershoot and the bottom curves give an example of 
both. 

Undershoot is evident is nearly all of the NPS-IRSTD images. The scan 
direction is from left to right across the scene, and the dark region to the right of 
the bright objects can be attributed to undershoot caused by the AC-coupling of the 
detector outputs. Comparing the prominence of the dark regions caused by AC- 
coupling effects to the influence of the optics and the detectors as demonstrated by 
the deconvolution routines points out the significant impact AC-coupling has on 


system operation. 


Be 


LOOHSYSONN ONY dOOYd 


LINDYID INITMNOI-JVY 
WOY¥d INdGLNO 





LINDYID INITdNAOI-JVY 


| | OL INdNl 


LOOHSYIONN 








LINDYIO INITdMAOI-OY 
WOH INdLNO 





{INDYIO YBNIMWNOBI-IJAV 


| Ol INdNl 


dood LiINDYIO BNITdNOI-OY 
WOHF ILNdINd 
INnHID ONITdAOI-OY 


1 
O1 I2NdNil 


Figure 46 AC-coupling Effects on Detector Output Signals 
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APPENDIX B 


IMAGING PRINCIPLES EXPRESSED THROUGH 
LINEAR SYSTEM THEORY 

This appendix is provided for the reader unfamiliar with the linear systems 
approach to image processing. The treatment of linear system theory here will be 
conceptual, and most general concepts will be simply stated without proof or 
argument. For a mathematical formulation of the topic the reader is referred to 
Lloyd [Ref. 9]. A mathematical treatment specifically directed toward the Fast 
Fourier Transform which was implemented in the image processing program can be 
found in Brigham [Ref. 14]. As in most situations where an assumption regarding 
the linearity of a process is made, there are circumstances in which non-linear 
effects will cause observations to deviate from predicted results. Lloyd [Ref. 9] 
points out the following situations in thermal imaging when linearity may be 
violated: 

* application of certain signal processing methods 

* given a spatially variant imaging process 

* when noisy conditions exist (system mapping is not single valued) 
The degree to which these non-linearities impact the imaging process varies widely 
between differing systems. 

A good starting point in the discussion of linear system theory is to introduce 
the concept of the Fourier function pair. Any periodic single valued function can 
be equivalently expressed as a weighted sum of a set of sine and/or cosine 
functions. A common example is a square wave, defined as a function of time and 
having a period T, which can be expressed as 


2unt 
qf 





eri 
Square Wave = > — cos ( ) (for n odd) . 
n 


a=} 
In this form each term in the sum represents the magnitude of a cosine wave with a 


frequency @=2nn/T. 
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The functions f(@) and f(t) from a unique pair. When f(t) is non-periodic, this 


relationship can often be expressed with Fourier integrals such that 


l fe —1twt 
ome [(0e dt 
PAE es 


f= | Fae! da, 


In shorthand notation ¥ [f(t)] defines the Fourier transform of f(t), then the Fourier 


pair can be expressed 

F(w)=F [f(t)] 

f(t)=F [F(o)] 
This is important because filters affect the various frequency components of a signal 
differently. If g(@) contains the information on the filter’s affect on the frequency 
components of the signal, then the output from the filter is 

h(w)=f(@) + g(@) 

In spatial coordinates f(x) and f(k) form the Fourier pair where k is related to 
spatial frequency as @ is to temporal frequency. 

In imaging systems the filters are the various parts of the system through 
which the signal must pass before it is seen by the observer; lenses, mirrors, 
obstructions, windows, detectors, electronics, and the display. One strength of the 
linear system approach is that each component can be treated as an independent 
filter. The final image is obtained by multiplying together the effects of all the 
components. 

Most of the fundamental ideas imbedded in linear filter theory can be 
illustrated using the imaging process for a simple lens as an example. The 
treatment for this discussion was taken from unpublished course notes written by A. 
W. Cooper and E. C. Crittenden, Jr. [Ref. 15]. Consider the scene to consist of a 
single point source located on the z-axis an infinite distance to the left of the x-y 
plane in Figure 47. The radiation from this source will arrive as plane waves 
parallel to the x-y plane. If the x-y plane contained a perfect lens of infinite extent, 


the light would be focussed to a point in the X-Y plane. 
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Point Spread Function 
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} Point Source 


Figure 47 Imaging of a Point Source at Infinity with a Circular Lens 
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However, since all lenses are of finite size, a point source will not be 
focussed to a point. Instead, the light will be distributed in a manner similar to that 
shown in Figure 48, which was taken from Hecht (Ref. 16]. The distribution 
represented by the curve in Figure 48 is referred to as the point spread function. 

The point spread function can be analytically determined by integrating the 
electric field incident on the x-y plane to find the resulting diffraction pattern in the 


X-Y plane, 


° 


F(X yar | | A (x,y) pik IRM Xx + Yy) dxdy 


In this expression A(x,y) contains the electric field amplitude and phase information 
at the lens, k=2n/A where A reperesents the wavelength of the light, R is the 
distance from the x-y plane origin, and E(X,Y) ts the electric field amplitude at the 
image plane. This expression can be put in the form of a standard Fourier integral 
by defining the radian spatial frequencies K,=(kX/R)=(2mX/AR) and 
(K,=kY/R=2nY/AR) so that 


i(K yx+K ) 


a y 
BK, Ky) = | | Atne Y" dxdy 


It can be seen that E and A are a Fourier transform pair. For the case of a point 
source on the optic axis an infinite distance from the optical system, the amplitude 
and phase function A(x,y) can be treated as a constant. Thus the distribution in the 
X-Y plane for a point source becomes just the Fourier transform of the shape of the 
optical component’s aperture. Figure 49, from Lloyd [Ref. 9], shows the two 
dimensional Fourier transform for several simple geometric shapes. 

A system for which the form of the point-spread function is the same for all 
points in the image is called spatially invariant. A linear system is described by 
behavior such that the scene can be decomposed into a sum of all its components. 
Therefore, for a linear, spatially invariant system the image projected onto the image 


plane is the result of convolving the point-spread function with each point in the 


scene. Using shorthand notation, a convolution is defined by the symbol so that 


h(x)=f(x)4 g(x) 
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f(x,y) 


F(k,) 





Figure 48 Point Spread Function for a Circular Aperture 
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GRAPHICAL REPRESENTATION FOURIER TRANSFORM 
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Figure 49 Two Dimensional Fourier Pairs 
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While computation of a convolution to evaluate the imaging characteristics of 
a system seems somewhat imposing, linear system analysis provides the mechanism 
for simplifying the problem. A fundamental theorem states that the Fourier 
transform of the convolution of two functions can be expressed as the product of 
the Fourier transforms of the functions. 

F ThOQ]=F [f@) xg l=F [f(x)] - F [g(x)]=F(k) » Gk) 
Application of this theorem involves moving the problem at hand back to the x-y 
plane where the functions are expressed in terms of frequencies. In other words, 
the convolution can be computed from the product of the Fourier transform of the 
scene distribution with the Fourier transform of the point spread function (which is 
the lens aperture shape). 

Extending the concept from one filter to many involves using the output from 
the last filter as the input to the next filter. Computation of the effect of the 
current filter is most easily accomplished in the frequency domain. Therefore, 
leaving the results from the previous filter in terms of frequency allows the effect of 
both filters to be computed as a product of the two. It can now be seen that 
removing the effect of a filter can be accomplished by dividing by its function in 
the frequency domain. This is known as deconvolution, and it is the process 
applied to remove the systems influence from image data. 

In practice the spatial fluctuations of scene brightness are described as a 
modulation. Modulation is sometimes defined as 

lia Ome tere) 
for both M,, the modulation of the image and M,, the modulation of the object, 
where I is the irradiance. From this it is straightforward to show that for a circular 
lens, higher spatial frequencies will have a much smaller modulation in the image 
than in the scene. This effect will be less noticeable for lower spatial frequencies. 


Figure 50 diagrams this concept. 
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Figure 50 Modulation of Point Source Distributions at Infinty Imaged onto the 
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The relationship of spatial frequencies and their corresponding degree of 
modulation is expressed in functional form by the modulation transfer function 
(MTF), where 

MTF = M,/M, 
Thus it is the MTF which is used to define the effect a system component will 
have on a signal passing through it. 

Application of Fourier optics theory provides a vehicle for showing that the 
MTF of an optical system can be theoretically derived to be the autocorrelation of 
the pupil function (described by the physical shape of the optical component) for a 
system imaging non-coherent light. This derivation is contained in Cooper and 
Crittenden’s Ph3252 course notes (Ref. 15]. Parker [Ref. 2] used ray tracing 
software as an analytical method to calculate the MTF of the NPS-IRSTD’s optical 
system, for both modified and unmodified versions. The results of these 


calculations [Ref. 2] are presented as Figure 51. 
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Figure 51 MTFs for the ADM and the NPS-IRSTD. (Top) AN/SAR-8 ADM 
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APPENDIX C 


I. COMPUTER PROGRAM LISTINGS 
BEGIN.FOR 


PROGRAM BEGIN.FOR WRITTEN BY: R. C. ENGEL 
LAST REVISED: 13 AUGUST 1989 
BY. C. ENGEL 


INTEGER*2 IBACK(15360), NUM(2), NUM1 (3) 

REAL XRES (90) 

DIMENSION IRED(0:255),RED(0:255), IGREEN (0:255) 
DIMENSION IBLUE(0:255),BLUE(0:255), GREEN (0:255) 


FORMAT STATEMENTS 
2 FORMAT (/) 


THESE ARE DT-IRIS CALLS TO INITIALIZE THE DRIVERS AND TO 
DISPLAY THE CONTENTS OF FRAME BUFFER O ON THE SCREEN. 


CALL ISINIT 
CALL ISOTFR(0) 
Can ISDISP (1) 


MENU SECTION. 


30 WRITE (*,2) 
WRITE (*,*) ‘1=CREATE COLOR SCALES FOR FRAMES 0 AND 1’ 
WRITE (*,*) ‘2=CHANGE OUPUT LOOKUP TABLES 2-7 
* (GAUSSIAN) ' 
WRITE (*,*) '3=CHANGE OUPUT LOOKUP TABLES 2-7 

* (LINEAR) ' 

WRITE (*,*) ‘4=CLEAR A DESIGNATED BUFFER’ 
WRITE (*,*) ‘5=VIEW A DESIGNATED BUFFER’ 
WRITE (*,*) '6=EXIT THE PROGRAM’ 
READ (*,*) NOPS 
GotN@PcesiT. 1 .OR. NCES w7Gr.)6)) Go TO 30 
Goro (99, 193, 201, 399,449,500) ,NOPS 


105 


ic 
LS. 
Cc 


THIS SECTION ADDS LABELS AND A COLOR SCALE ACROSS THE 
BOTTOM OF THE DISPLAY SCREEN FOR BUFFERS QO AND 1. 


99 DO 100 IFR=0,1 
DO 110 I=0,29 
DO 120 J=0,510,2 
IBACK (I*512+J+1) =J3/2 
IBACK (1*512+J+2) =J3/2 


120 
16 


CONTINUE 
CONTINUE 
CALL ISPUTP 
CALL ISSFNT 
NUM (1) =50 
NUM (2) =53 
CALL ISGPOS 
CALL ISTEXT 
NUM (1) =53 
NUM (2) =48 
CALL TSGPCS 
CALL ISnhocr 
NUM (1)=55 
NUM (2) =53 
CALL ISGPOS 
CALL ISTEXT 


NUM1 (1) =49 
NUM1 (2) =48 
NUM1 (3) =48 


CALL ISGPOCS 
CALL ISTEXT 
NUM1 (1) =49 
NUM1 (2) =50 
NUM1 (3) =53 
CALL ISGPOS 
CAGie lols... 
NUM1 (1) =49 
NUM1 (2) =53 
NUM1 (3) =48 
CALL ISGPOS 
CALL ISTEXT 
NUM1 (1) =49 
NUM1 (2) =55 
NUM1 (3) =53 
CALL ISGPOS 
CALL ISTEXT 
NUM1 (1) =50 
NUM1 (2) =48 


(tPA S0, 0, 15300, TBACK) 
(0) 


(460, 45) 
(IFR, 2, NUM) 


(460, 95) 
(IFR, 2, NUM) 


(460,145) 
(IFR, 2, NUM) 


(460,192) 
(IFR, 3,NUM1) 


(460, 242) 
(IFR,3,NUM1) 


(460, 292) 
(IFR, 3,NUM1) 


(460, 342) 
(IFR,3,NUMi) 
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C 
C 
C 


c 
c 
Se 


Cc 
C 


Cc 


NUM1 (3) =48 
CALL ISGPOS 
CALL ISTEXT 
NUM1 (1) =50 
NUM1 (2) =50 
NUM1 (3) =53 
CALL ISGPOS 


(460,392) 
(Cla 2, NUM) 


(460, 442) 


201 WRITE 


CALL ISTEXT 
100 CONTINUE 


GOVTO, 20 


THIS SECTION BUILDS OUTPUT TABLES 


(IFR,3,NUM1) 


(GAUSSIAN) 


NUMBER OF THE OUTPUT TABLE TO 


’2-5 INCLUSIVE ALLOWED’ 


198 WRITE (*,2) 
WRITE (*,*) ‘ENTER THE 
* MODIFY’ 
WRITE (*,*) 
READ (*,*) KTAB 
WRITE (*,*) ‘ENTER THE 
READ (*,*) RMEAN 
WRITE (*,*) "ENTER THE 
READ (*,*) RDEV 
WRITE (*,*) "ENTER THE 
READ (*,*) GMEAN 
WRITE (*,*) ‘ENTER THE 
READ (*,*) GDEV 
WRITE (*,*) ‘ENTER THE 
READ (*,*) BMEAN 
WRITE (*,*) ‘ENTER THE 
READ (*,*)  BDEV 


CALCULATE THE VALUES FOR 


BO 199 T—0, 255 


IRD=INT (255* (EXP (-0. 
IGR=INT (255* (EXP (-O. 
EPBE=INT(255* (EXP (-0. 


MEAN VALUE FOR RED’ 

STD DEVIATION FOR THE RED’ 
MEAN VALUE FOR GREEN’ 

STD DEVIATION FOR THE GREEN, 
MEAN VALUE FOR BLUE’ 


STD DEVIATION FOR THE BLUE’ 
THE TABLE 


BG caren ies 2 (RDEV) **2))) ) 
5* ( (I-GMEAN) **2) / ( (GDEV) **2)))) 
2 ( (1 BMBAN) <*2)/ (@BDEV) **2)))) 


CALL) ISLDOV (KTAB, 1, IRD, IGR, 1BL) 


199 CONTINUE 


CALL ISOUTS (KTAB) 


COLT 30 


READ 


iar) 


(25%) 


‘WHICH TABLE DO YOU WANT TO CHANGE 


NTAB 


BUILD A COLOR TABLE USING THREE POINTS TO DEFINE TWO 
LINES REPRESENTING THE VALUE FOR THE GIVEN INDEX. 


22 oO) 
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IF (NTAB .LI. 2° 7OR. NPAB Gl ope Gcoetond 
WRITE (*,*) "ENTRIES MUSTOALL BE INTEGERS EROMeO TO 


* 255! 
e 
C CALCULATE THE RED LUT VALUES 
G 
200 WRITE (*,*) ’ENTER THE LOW RED VALUE INDEX, 
x MAGNITUDE’ 
READ (*,*) LWRDIDX, LWRDVAL 
IF (LWRDIDX .GT. 255 .CR. LWRDVAL) .Gi. 255) Gomro 200 
205 WRITE (*,*) ‘ENTER THE MIDDLE RED VALUE INDEX, 
* MAGNITUDE’ 
READ (*,*) MDRDIDX, MDRDVAL 
IF (MDRDIDX .GT. 255 .OR. MDRDVAL [GT. 255) Gom@o 206 
IF (MDRDIDX .LE. LWRDIDX) THEN 
WRITE (*,*) “INDEX TOO LOW’ 
GO TO 205 
ENDIF 
210 WRITE (*,*) ‘ENTER THE HIGH RED VALUE INDEX, 
* MAGNITUDE’ 
READ (*,*) IHRDIDX, IHRDVAL 
IF (IHRDIDX .GT. 255 .OR. IHRDVAL .GT. 255) GO TO 200 
IF (IHRDIDX .LT. MDRDIDX) THEN 
WRITE (*,*) ‘INDEX TOO LOW’ 
GS WO 210 
ENDIF 
IRED(LWRDIDX) = LWRDVAL 
IRED (MDRDIDX) = MDRDVAL 
IRED(IHRDIDX) = IHRDVAL 
IF (LWRDIDH -5O. ©) Ge: To 216 
DO 215 I=0, LWRDIDX-1 
IRED (I) =0 
215 CONTINUE 
216 SLP=1.0* (MDRDVAL-LWRDVAL) / (MDRDIDX-LWRDIDX) 
DO 220 I=LWRDIDX+1, (MDRDIDX-1) 
RED (I) =LWRDVAL+SLP* (I-LWRDIDX) 
IRED (I) =INT (RED (I) ) 
220 CONTINUE 
IF (IHRDIDX .EQ. MDRDIDX) GO TO 230 
SLP=1.0* (IHRDVAL-MDRDVAL) / (IHRDIDX-MDRDIDX) 
DO 225 I=MDRDIDX+1, (IHRDIDX-1) 
RED (I) =MDRDVAL+SLP* (I-MDRDIDX) 
IRED (I) =INT (RED (1) ) 
225 CONTINUE 
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€ CALCULATE THE GREEN LUT VALUES 


G 
230 


AS) 


240 


245 
246 


250 


2a0 


260 


WRITE (*,*) ‘ENTER THE LOW GREEN VALUE INDEX, 
MAGNITUDE’ 

READ (*,*) LWGRIDX, LWGRVAL 

IF (LWGRIDX .GT. 255 .OR. LWGRVAL .GT. 255) GO TO 230 

WRITE (*,*) ‘ENTER THE MIDDLE GREEN VALUE INDEX, 
MAGNITUDE’ 

READ (*,*) MDGRIDX, MDGRVAL 


IP SMOGRIDX |G 255 .OR. MDGRVAL).GT: 255) Go TO 235 
IF (MDGRIDX .LE. LWGRIDX) THEN 
WEEE (+, ~)) INDEX TOO Low? 


GOero 255 
ENDIF 
WREEE(*,“~) “ENTER THE HIGH GREEN VALUE INDEX, 
MAGNITUDE’ 
READ (7) = ) IHGRIDX, IHGRVAL 


DESGLHGR ED MGT. 255 .OR. TBGRVAL .GT. 255) GO TO 240 
tP SCIHGRIDX .LEo MUGRIDX) THEN 
WRITE (*,*) ‘INDEX TOO LOW’ 


GO TO 240 
ENDIF 
IGREEN (LWGRIDX) = LWGRVAL 
IGREEN (MDGRIDX) = MDGRVAL 


IGREEN (IHGRIDX) = IHGRVAL 
fh WweRmDx FO. 0) GO TO 246 
DO 245 I=0, LWGRIDX-1 
IGREEN (I) =0 
CONTINUE 
SLP=1.0* (MDGRVAL-LWGRVAL) / (MDGRIDX-LWGRIDX) 
DO 250 I=LWGRIDX+1, (MDGRIDX-1) 
GREEN (I) =LWGRVAL+SLP* (I-LWGRIDX) 
IGREEN (I) =INT (GREEN (TI) ) 
CONTINUE 
IF (THGRIDX (EO. MDGRIDX)) GO TO 260 
SLP=1.0* (IHGRVAL-MDGRVAL) / (IHGRIDX-MDGRIDX) 
DO 255 I=MDGRIDX+1, IHGRIDX-1 
GREEN (1) =MDGRVAL+SLP* (I-MDGRIDX) 
IGREEN (I) =INT (GREEN (I) ) 
CONTINUE 


CALCULATE THE BLUE LUT VALUES 


WRITE (*,*) ‘ENTER THE LOW BLUE VALUE INDEX, 
MAGNITUDE’ 
READ (*,*) LWBLIDX, LWBLVAL 
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aqaaa 


IF (LWBLIDX .GT. 255°OR. LWBLVAM@e GT, 255) 9"Goero 260 


265 WRITE (*,*) ‘ENTER THE MIDDLE BLUE VALUE INDEX, 
* MAGNITUDE’ 
READ (*,*) MDBLIDX, MDBLVAL 


IF (MDBLIDX .GT. 255 .OR. MDBLVAL .GT. 255) GO TO 265 
IF (MDBLIDX .LE. LWBLIDX) THEN 
WRITE (*,*) ‘INDEX TOO LOW’ 


GO TO 265 
ENDIF 
270 WRITE (*,*) ’ENTER THE HIGH BLUE VALUE INDEX, 
* MAGNITUDE ’ 
READ (*, 7) IHBLIDX, IHBLVAL 


IF (IHBLIDX .GT. 255 OR. IHBLVAL [Gi 255) Goel 0 270 
IF (L2HBLIDX (Li MeBLID> Mirae 
WRITE (*,*) ‘INDEX TOO LOW’ 


GOvlGr27 2 
ENDIF 
IBLUE (LWBLIDX) = LWBLVAL 
IBLUE (MDBLIDX) = MDBLVAL 
IBLUE(IHBLIDX) = IHBLVAL 


IF (LWBLIDX GEO. 0) Couto cic 
DO 275 I=0, LWBLIDX-1 
IBLUE (I) =0 
275 CONTINUE 
276 SLP=1.0* (MDBLVAL-LWBLVAL) / (MDBLIDX-LWBLIDX) 
DO 280 I=LWBLIDX+1, (MDBLIDX-1) 
BLUE (I) =LWBLVAL+ (SLP* (I-LWBLIDX) ) 
IBLUE (I) =INT (BLUE (I) ) 
280 CONTINUE 
IF (IHBLIDX .EQ. MDBLIDX) GO TO 290 
SLP=1.0* (IHBLVAL-MDBLVAL) / (IHBLIDX-MDBLIDX) 
DO 285 I=MDBLIDX+1, IHBLIDX-1 
BLUE (I) =MDBLVAL+SLP* (I-MDBLIDX) 
IBLUE (I) =INT (BLUE (I) ) 
285 CONTINUE 


DT-IRIS CALLS TO PUT THE CALCULATEDS VAEUES IN TTHE TAELE. 


290 CALL ISDISP (1) 
DO 296 I=0,255 
CALL ISLDOV(NTAB,I, IRED (I), IGREEN(I) , IBLUE (I) ) 
296 CONTINUE 
CALL ISOUTS (NTAB) 
GO TO 30 
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THIS SECTION CLEARS A DESIGNATED FRAME BUFFER. 


399 WRITE (*,*) ‘WHICH FRAME BUFFER WOULD YOU LIKE TO 


i CLEAR?’ 
aaVANID) (25, BS) MF RAM 
CALL ISFCLR (MFRAM) 
GOTO 30 


THIS SECTION DISPLAYS A DESIGNATED FRAME BUFFER. 


449 WRITE (*,*) ‘WHICH FRAME BUFFER WOULD YOU LIKE TO 


x VIEW?! 
READ (*,*) MFRAM 
CALL ISOTFR (MFRAM) 
GOTO 30 


PROGRAM ENDS. ISEND MUST BE USED TO PREVENT POSSIBLE 
CONLICTS IN MEMORY. 


500 CALL ISEND 


STOP 
END 


DMA.FOR 


PROGRAM: DMA.FOR IPSILON ACS IR, (Co TANGINE, 
LAST REVISED: 07 AUGUST 1989 
BY: RAY ENGEL 


THIS PROGRAM IS FOR CONTROLLING DMA INPUT TO RAM FROM 
70000H TO AOOOO0H. DATA RECEIVED BY DMA TRANSFER IS THEN 
SORTED TO RAM FROM 37CO0OH TO 70000H BEFORE BEING WRITTEN 
IN BLOCKS TO FRAME BUFFER MEMORY AT AOOOOOH. 

INTEGER*4 ISADRZ, IDADR, ISADR, IRADD1 

INTEGER*4 IRADD2,IRADD3, IDADRX 

INTEGER*2 NADDR(90),GDT1 (24) 
FORMAT STATEMENTS. 
2 FORMAT (/) 
MENU SECTION. 
50 WRITE (*,2) 

WRITE (*,2) 
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= SINGLE DMA SET AFTER SECTOR FIND’ 
IMMEDIATE SINGLE DMA SET’ 
EXIT THE PROGRAM’ 


WRITE (*,*) ’ 
WRITE (*,*) ° 
WRITE (*,*) ° 
WRITE (*,2) 
WRITE (*,*) ' ENTER VALUE 1-3 ’ 

READ (*,*) IMODE 

IF (IMODE .LT. 1 JOR. IMOPEs- GT. Ss) GO 106850 
IF (IMODE .EQ. 2) GO TO 120 

IF (IMODE .EQ. 3) GO TO 500 


WNrR 
ll 


THIS SECTION FINDS THE END OF ROTATION PULSE AND THEN 
COUNTS 2.66 SECONDS PER SECTOR TO FIND THE DESIRED AREA 
OF OBSERVATION. THE FIRST STEP IS TO PUT PORTS A AND B 
IN INPUT MODE. SEE PAGES 22 AND 24 OF THE PDMA-16 
MANUAL [REF. 4]. 


SET PORT B FOR INPUT. 
CALL OPORTB (754, 0) 


INPUT FROM KEYBOARD THE SECTOR DESIRED TO DISPLAY 


60 WRITE (*,*) “ENTER THE SECTOR TO PESPEAY, (052-2) 


ALLOWED) ’ 

READ (*,*) SCrR 

IF (SCTR LT. 022 OR. SCIR Giles Jes) meco  LOmne 
TIMES 275607 5¢C0r 


NEXT THE END OF ROTATION SIGNAL MUST BE FOUND. THIS IS 
ON THE DO BIT OF PORT B AND IS MUCH LONGER THAN THE 
SECTOR SIGNAL, WHICH IS ALSO ON DO. WHEN DO GOES LOW, 
THE 30 KHZ ROTATION SIGNAL IS WATCHED FOR SIX CYCLES, 
THEN DO IS CHECKED AGAIN. IF DO IS STILL LOW, END-OF 
ROTATION HAS BEEN FOUND. 


70 CALL IPORTB(753,ICLKS) 
IEOR = IAND(ICLKS,1) 
IF (IEOR.EOQ.1) GO To 70 
DO 80 I=1,6 
81 CALL IPORTB(753, ICLKS) 
IEOR = IAND(ICLKS, 2) 
IF (IEOR.EQ.0) GO TO 81 
82 CALL IPORTB(753, ICLKS) 
IEOR = IAND(ICLKS,2) 
ff (1EOR.£O.2) GO TO 92 
80 CONTINUE 
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DO WENT LOW AND Dl HAS CYCLED SIX TIMES, IF DO IS STILL 
LOW THEN EOR HAS BEEN FOUND. 


CALE EPORTE (753, L@LKS) 

IEOR = IAND(ICLKS,1) 

TE MGhEOR. EO. 1) “GO TO 70 

CAULPGELEIM (IHR, IMIN, ISEC, LAN) 
CLKTIM=IHR*3600+IMIN*60+ISEC+IHN*0.01 

WE iia ~e  BORRE OUND LOOKING FOR SECTOR *, S€TR 


THE SECTOR IS FOUND BY COUNTING 2.66 SEC/SECTOR FROM THE 
PC’S SYSTEM CLOCK. THIS PROCEDURE WAS NECESSARY BECAUSE 
JITTER IN THE SECTOR CLOCK WOULD INADVERTENTLY TRIGGER 
CYCLES WHICH DISRUPT THE COUNTING PROCESS. 


90 CALL GETTIM(IHR, IMIN, ISEC, IHN) 


CLATIMP-THR* 3600+ IMIN“60+1SEC+IHN*0.01 
cc hLATiMi-ChLATIMy; “Gt. TIME) GO TO 120 
GO TO 90 


HAVING COMPLETED THE PREVIOUS LOOP, THE CORRECT SECTOR 
HAS BEEN LOCATED. THE NEXT THING TO DO IS FIND THE 
SlAtvsser THE 30 KHa SCAN CLOCK.” SINCE IT WILL HAVE 
JUST GONE HIGH AT THE END OF A DMA TRANSFER, I NEED TO 
ENTER THE DMA LOOP JUST AFTER A LOW TO HIGH TRANSITION. 


POR cALl, [PORTE (753, ICLKS) 


ISCAN = IAND(ICLKS, 4) 
PE ESeAN BO. 4) GO TO 120 


125 CALL IPORTB(753, ICLKS) 


ISCAN = IAND(ICLKS, 4) 
iP S(lSCANSEOQ.0) GO TO 125 


RHEE TesceerlON SETS UP AND EXECUTES 3 CONSECUTIVE DMA 
TRANSFERS FOR 65528 BYTES EACH TO STARTING ADDRESSES 
H7000° TO H9OO0O0. 


BO 150° 1=7,9 
CALL OPORTB(754,0) 


DISABLE ANY DMA OPERATION IN PROGRESS. 


CALL OPORTB (10,5) 


SET THE PAGE (64 KB REGION) FOR THE TRANSFER 
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CALL OPORTE (2a ais) 


LOAD TRANSFER OFFSET = 0, ALL TRANSFERS START ON 
A BOUNDARY. 


CALL OPORTB (2, 0) 
CALL OPORTB (2,0) 


LOAD THE BYTE COUNT REGISTERS, LOW Sr le then enews Te 
PDMA 16 ALWAYS DOES 1 EXTRA TRANSFER. 

NUMBER = 256*HB + LB + 1 WHERE HB AND LB ARE 2ND 
ARGUMENTS. 


CALL OPORTB (3,247) 
CALL OPORTB (3,255) 


SET TRANSFER DIRECTION AND MODE (PAGE 13 OF I 9) MANUAL). 
CALL OPORTB (11, 69) 


SCAN BIT SHOULD STILL BE HIGH, LOOK FOR It) tO .G ete: 
WHEN IT GOES HIGH AFTER THAT, START THE DMA. 


LS) CALL LPORTE(? Sal ebis) 


ISCAN = IAND(ICLKS, 4) 
ER (ESCAN EO 4a GOuTomi51 


152 (CALL LPOREB (7 S27 FCLNs) 


ISCAN = IAND(ICLKS, 4) 
IF (ISCAN.EQ.0) GO TO 152 


START THE DMA. BIT D7 OF DMA REG, SEE PAGE 24 OF MANUAL. 
CALL OPORTB (754,128) 
ENABLE 8237 MASK BIT 


CALL OPORTB(10,1) 
PRINT *,’DMA TRANSFER RUNNING’ 


THIS SECTION CHECKS FOR COMPLETION OF DMA TRANSFER 


155 CALEMIPORTE (3,11) 


CALL IPORTB (3, IH) 
NLEFT = 256*IH + IL 
IF (NLEFT.NE.65535) GO TO 155 


150 CONTINUE 
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THIS PROGRAM MOVES DATA FROM STANDARD RAM TO 
FRAMEGRABBER MEMORY FOR DISPLAY. 


THE DATA SEQUENCE IS THE PHYSICAL LOCATION OF THE NEXT 
CONSECUTIVE DETECTOR DATA IN RAM. (ie. DETECTOR 1 IS IN 
BASE + 0 WHILE DETECTOR 2 IS BASE + 6.) 

DATA NADDR/24, 36, 48,60, 72, 84, 96,108,120,132,144,156, 
plage 2, 25,31,49,61, 73,85, 97,109, 1217133/7145,157, 
Pio eels 2G, 30750, 62, 74,60, 96,110, 122,134,146, 158, 
Ome?) 39, 51,63, 75,87, 99, 111,123,135, 147,159, 
iio, 5,260,400, 02,64) 76,88, 100, 112, 124,136, 148,160, 
ii 2eoelG, 29, 41,5365, 77,89, \0leils, 125,037,149, 161, 
iG BG Gio allet 78 


SCREEN MESSAGE. 


WRITE (*,2) 

WRITE (*,*) ‘ DMA TRANSFERS COMPLETE. UNSCRAMBLING 
* DATA AND ’ 

WRITE (*,*) ’ WRITING TO FRAMEGRABBER MEMORY ’ 
WRITE (*,2) 


UNSCRAMBLING FROM DMA TRANSFERS TO 70000H AND 80000H FOR 
IMAGE DESTINED FOR FRAME BUFFER 0. 

ISADR* IS REFERRING TO SOURCE ADDRESSES. 

IDADR* REFERS TO DESTINATION ADDRESSES. 


ISADRZ=458752 
I=0 
205 CALL PEEK(IOFST, ISADRZ+I,1) 
I=I+1 
Pe aChOnsT oNE.0) GO TO 205 
ISADRZ=458757+I1 
IDADR=228352 
BO 210 1=1, 90 
N=NADDR (I) 
DO 220) 0=0, 363 
ISADR=ISADRZ+ (180*J+N) 
DO 230 L=0,4 
IDADRX=IDADR+L*512 
CALL BLKMOV (ISADR, IDADRX, 1) 


250 CONTINUE 
IDADR=IDADR+1 
220 CONTINUE 


IDADR=IDADR+2196 
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210 CONTINUE 
ISADRZ=524288 
I=0 
245 CALL PEEK(IOFST, ISADRZ+I,1) 
I=I+1 
IF (1OFST.NE.O0) GO TO 245 
ISADRZ=524293+I 
IDADR=228352+364 
DO 250 I=1,90 
N=NADDR (I) 
pe) CAS) ate 7 
ISADR=ISADRZ+(180*J+N) 
DO 270 L=0,4 
IDADRX=IDADR+L*512 
CALL BLKMOV (ISADR, IDADRX,1) 


270 CONTINUE 
IDADR=IDADRt1 
260 CONTINUE 


IDADR=IDADR+2412 
250 CONTINUE 


BLOCK TRANSFERS TO FRAME BUFFER 0. 


DO 290 I=0,3 
IRADD=10485760+1*57600 
ISADR=228352+1*57600 
IER = BLKMAT(GDT,ISADR, IRADD, 28800) 
290 CONTINUE 


SORTING FROM 2ND AND 3RD DMA TRANSFERS FOR IMAGE 
DESTINED TO FRAME BUFFER 1. 


ISADRZ=524288 
I=0 
305 CALL PEEK(IOFST, ISADRZ+I,1) 
I=I+1 
IF (IOFST.NE.0) GO TO 305 
ISADRZ=524293+I1 
IDADR=228352 
DO 310 I=1,90 
N=NADDR (I) 
DO 320 J=148,363 
ISADR=ISADRZ+ (180*J+N) 
DO 330 L=0,4 
IDADRX=IDADR+L*512 
CALL BLKMOV (ISADR, IDADRX, 1) 
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330 CONTINUE 
IDADR=IDADR+1 
320 CONTINUE 
IDADR=IDADR+2344 
310 CONTINUE 
ISADRZ=589824 
I=0 
S45 GAml, PEEK {(TOPST, ISADRZ+1-1) 
I=I+1 
IF (IOFST.NE.0) GO TO 345 
ISADRZ=589829+I 
IDADR=228352+216 
Dows5o) 11,90 
N=NADDR (I) 
DO 360 J=0,295 
I SADR=ISADRZ+ (180*J+N) 
DO 370 L=0,4 
IDADRX=IDADR+L*512 
CALL BLKMOV (ISADR, IDADRX, 1) 


370 CONTINUE 
IDADR=IDADR+1 
360 CONTINUE 


IDADR=IDADR+2264 
350 CONTINUE 


BLOCK TRANSFERS TO FRAME BUFFER 1. 


DO 3005-0, 3 
IRADD=10485760+1*57600+262144 
ISADR=228352+1I*57600 
IER = BLKMAT(GDT, ISADR, IRADD, 28800) 

390 CONTINUE 


COyTO: 50 
PROGRAM ENDS. 
500 3STOP 
END 
IMAGE.FOR 
PROGRAM: IMAGE.FOR WRETIENe Sys. C. ENGEL 
LAST REVISED: 18 AUGUST 1989 
Bye Ro CC. ENGEL 
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IMPLICIT INTEGER*2 (I-N) 


IMPLICIT COMPLEX (C) 
DIMENSION CG(31zZ) , REEINE (912)), REI(512)9 XRES (20) 
DIMENSION AMAG (512) ,S8 (512) 
INTEGER*2 LTROW1 (0:511) , TROW2 (0: S511)> TROWS (Oe 1) 
INTEGER* 2 TROW4 (02511), LROWS (03511), LINE (32) 
INTEGER* 2 IBACK (15360) ,NUM(2) , NUM1 (3) , IBTTM(9) 
INTEGER* 2 ISIXA(0:255),ISIXB(0:255) , NFRAME, ITOP (9) 
CHARACTER*30 FILES, ANS$*1 
€ 
Cc FORMAT STATEMENTS. 
c 
2 FORMAT (/) 
i 
Cc ISINIT MUST BEGIN ALL DT-IRIS PROGRAMS. Devise euUie 
Cc THE FRAMEGRABBER IN DISPLAY MODE. 
G 
CALL ISINIT 
CALL ISbiISP (1) 
ic 
Cc MENU SECTION. 
Cc 
50 WRITE (*,2) 
WRITE (*,*) ‘ 1 = select the frame buffer to display ’ 
WRITE (*,*) ’ 2 = select/change output look-up table ’ 
WRITE (*,*) ’ 3 = make 6 degree image in buf 2 ’ 
WRITE (*,*) * 4 = correct for detector responsi ty, 
“ differences ’ 
WRITE (*,*) ’ 5 = fill in lines for bad detectors ’ 
WRITE (*,*) 0’ 6 = apply MIPWcorrection,” 
WRITE (*,*) ' 7 = smooth out 1 detector to 5 row 
2S expansion 
WRITE (*,*) ' 8 = store a frame to a file ’ 
WRITE (*,*) ' 9 = retrieve a file to a buffer ’ 
WRITE (*,*) ‘10 = copy one buffer to another ’ 
WRITE (*,*) ’11 = subtract one buffer from another ’ 
WRITE (*,*) ‘12 = 6xit this program ~ 
READ (*, *) KSELCT 
IF AKSELCT .LI. 1 [,OR-SKSELCT Gi eke) = COmLe moe 
GO TO (100, 200, 300, 400,500, G00, 7007800, 900, 
*V000;7 11007 0500), KSELET 
€ 


Cc THIS SECTION SELECTS THE FRAME BUFFER TO DISPLAY 
c 


100 WRITE (*,*) ’WHICH FRAME BUFFER WOULD YOU LIKE TO 
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. VIEW?’ 
READ (*,’ (I2)’) NFRAME 
CALL ISOTFR (NFRAME) 
GO TO 50 


THIS SECTION SELECTS THE ACTIVE OUTPUT LOOK-UP TABLE. 


200 WRITE (*,*) ' SELECT LUT BY NUMBER 0-5’ 
READ (*,*) ITAB 
CALL ISOUTS (ITAB) 
GO TO 50 


THIS SECTION BUILDS A 6 DEGREE HORIZONTAL IMAGE IN FRAME 


BUFFER 2 FROM THE 3 DEGREE IMAGES IN BUFFERS 0 AND 1. 


BOO CALL ESDIVC (0,2, 2) 
CAlnn LSDeve (1, 273) 
DO 310 J=0,449,5 
GALL ISGETP(2,J,0,5U2z,2ROW!) 
GAlu ISGHIP(3,3,0,512, LROWZ) 
BOes 20) K=0,511,2 
ISIXA (K/2) =IROW1 (K) t+IROW1 (K+1) 
ISIXB(K/2) =IROW2 (K) +IROW2 (K+1) 
320 CONTINUE 
DO 321 K=0,4 
L=J+K 
Chih © SPUTP(2, L, 0,256, 1S1XA) 
Bali ISPUTP(2,L, 256,256, 1S1XB) 
321 CONTINUE 
310 CONTINUE 


THIS SECTION ADDS LABELS AND A COLOR SCALE ACROSS THE 


BOTTOM OF THE DISPLAY SCREEN FOR FRAME BUFFER 2. 


NFRAME=2 
DO 330 I=0,29 
DO 340 J=0,510,2 
IBACK (1*512+J+1) =J/2 
IBACK (1*512+3+2) =J/2 
340 CONTINUE 
330 CONTINUE 
CALL ISPUTP (NFRAME, 450,0,15360, IBACK) 
CALL ISSFNT (0) 
NUM (1) =50 
NUM (2) =53 
CALL ISGPOS (460, 45) 
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CALL ISTEXT (NFRAME, 2, NUM) 
NUM (1) =53 

NUM (2) =48 

CALL ISGPOS (460,95) 

CALL ISTEXT (NFRAME, 2, NUM) 
NUM (1) =55 

NUM (2) =53 

CALL ISGPOS (460,145) 

CALL ISTEXT (NFRAME, 2, NUM) 
NUM1 (1) =49 

NUM1 (2) =48 

NUM1 (3) =48 

CALL ISGPOS (460,192) 

CALL ISTEXT (NFRAME, 3,NUM1) 
NUMI1 (1) =49 

NUM1 (2) =50 

NUM1 (3) =53 

CALL ISGPOS (460,242) 

CALL ISTEXT (NFRAME, 3,NUM1) 
NUM1 (1) =49 

NUM1 (2) =53 

NUM1 (3) =48 

CALL ISGPOS (460,292) 

CALL ISTEXT (NFRAME, 3, NUM1) 
NUM1 (1) =49 

NUM1 (2) =55 

NUM1 (3) =53 

CALL ISGPOS (460,342) 

CALL ISTEXT (NFRAME, 3, NUM1) 
NUM1 (1) =50 

NUM1 (2) =48 

NUM1 (3) =48 

CALL ISGPOS (460,392) 

CALL ISTEXT (NFRAME,3,NUM1) 


NUM1 (1) =50 
NUM1 (2) =50 
NUM1 (3) =53 


CALL ISGPOS (460,442) 

CALL ISTEXT (NFRAME, 3, NUM1) 
CALL ISOTFR(NFRAME) 

GO TO 50 


THIS SECTION CORRECTS FOR THE DIFFERENT RESPONSES OF THE 
DETECTORS IN THE ARRAY. 


THIS DATA SET CONTAINS THE VALUES OF RESPONSIVITY FOR 
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C THE LEAD ARRAY AS MEASURED BY AYERS IN HIS EVALUATION OF 

C THE NPS-IRSTD. 

eS 
BAA ARES 152271350, lag580,1.904, 1.900,1.747, 1.681, 
ZO 70g Colin g eer diesiG? 1° 631.583, 1-790, 0.400. 
aieerolewese, |. 782,1.702,1.691,1.750,1.763,1.612, 
Tle oi 1643, 1).592,1,642,1,661,1-581.0.,050)17611) 
mlleoo2, 1. 702, 1.722,1,800, 1,830, 1.650,1.381, 1.680, 
el G04), 7S, 1756-1) 7121). 693,1.431,, 1. ado, 
Pieise ot 0.060, 1uo65, 146991 794,1.835, 1.79% 
pete le 129) 1s 70), e791) 7297 1. 690, 1.679, 1.673, 
Pico. 142 120091668, 1.733), 1.820,1.720, 1.770% 
Selle 824, 9833-05035, 0.030,1.829, 0.030, 1.902, 
eo Opeee 7 G7 WC S4 te 04, 1.62171.711,1.751,1.662, 
*1.662,1.640,1.681/ 

@ 

C SELECT FRAME TO CORRECT. 

c 

400 WRITE (*,*) ' ENTER NUMBER OF FRAME TO CORRECT ’ 

READ (*,’(I2)’) NFRAME 

e 

C COMPUTE THE AVERAGE OF THE RESPONSIVITY VALUES. 

c 


AVGRES=0 
DO 410 I=1,90 

IF (XRES(I) .LT. ..7*AVGRES/I) XRES(I)=AVGRES/ (I-1) 
* AVGRES=AVGRES+XRES (I) 


410 CONTINUE 
AVGRES=AVGRES/ 90 


CORRECT EACH ROW OF VALUES BASED ON THE RELATIVE 
REE SONSIVITY., 


angaaAA 


DO 420 I=0,445,5 
CALL ISGETP (NFRAME,1I,0,512, I1ROW1) 
DO 430 J=0,511 
ITROW2 (J) =INT ( (IROW1 (J) -128) *AVGRES/ 
* MRES( CLASH Pl) ti26) 
IF (IROW2(J) .GT. 255) IROW2(J)=255 
430 CONTINUE 
DO 440 K=0,4 
CALL ISPUTP (NFRAME, I+K,0,512, IROW2) 
440 CONTINUE 
420 CONTINUE 
CALL ISOTFR(NFRAME) 
CALL 1SDISh ay 


121 


CLOVOCOrOter@ 


Qqaagaa GEOG) 


ORO@nG 


AAO eea © 


GOVTO 50 
THIS PROGRAM FILLS THE ROWS FROM BAD DETECTORS IN THE 
SCENE BY A LINEAR APPROXIMATION BETWEEN THE VALUES OF 
THE ROWS ABOVE AND BELOW THE BAD DETECTOR. 
ROWS BOUNDING BAD DETECTORS. 


DATA ITOP/44,59,144,199,219,294, 369, 404, 439/ 
DATA IBTTM/50, 80,150, 205,225, 305, 380,410, 445/ 


SELECT THE FRAME TO SE) CORRECTED: 


500 WRITE (*,*) ‘ ENTER THE NUMBER OF THE FRAME FILL IN ’ 


READ (*,’(I2)’) NFRAME 


DISPLAY THE FRAME BEING CHANGED. 
CALE TSBESECa) 
CALL ISOTFR (NFRAME) 
CALCULATE THE APPROPRIATE PIXEL VALUES AND WRITE THEM TO 
MEMORY. 
DO 505 I=1,9 
CALL ISGETP (NFRAME, ITOP (I),0,512, IROW1) 
CALL ISGETP (NFRAME, IBTTM(I),0,512, IROW3) 
NDIF=IBTTM (I) -ITOP (I) 
DO 510 J=ITOP (I) +1, IBTTM(I)-1 
DO 515 K=0,511 
IROW2 (K) =INT (1.0* ( (IROW3 (K) -IROW1 (K) ) /NDIF) * 
(J-ITOP (I) )) +IROW1 (K) 
515 CONTINUE 
CALL ISPUTP (NFRAME,J,0,512, IROW2) 
510 CONTINUE 


505 CONTINUE 


GoUrloysy 


THIS SECTION ALLOWS CORRECTIONS FOR THE DEGRADATIONS IN 
THE IMAGE CAUSED BY THE MTF OF THE LENS AN ae x ANE 2 
DIRECTION, THE CONVOLUTION CF These e ie eCrOnRe sean Thy THe en 
DIRECTION, AND THE DETECHGCR Mir IN @tHe YevrRaecrion. 


SELECT foe DESIRED BRAME TO CORRECT. 


600 WRITE (*,*) ‘WHICH FRAME TO MIF CORRECT ” 
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READ (*,' (12)"’) NERAME 

EE ONERAME -LT. 0 -OR. NFRAME .GT. 15) GO TO 600 
CALL ISDISP (1) 

CALL ISOTFR (NFRAME) 


c 
C CHOOSE THE PARTICULAR ASPECT TO MANIPULATE. 
C 
601 WRITE (*,*) ’ 1 = X-DIRECTION SCAN DECONVOLUTION ' 
WRITE (*,*) ‘ 2 = X-DIRECTION OPTICS ’ 
WRITE (*,*) ‘ 3 = Y-DIRECTION OPTICS ’ 
WRITE (*,*) ‘ 4 = RETURN TO MAIN MENU ° 
READ (*,*) ISELCT 
ty isk wene Url) .OR, ISELCT -Gr. 5) GO TO.605 
GO TO (602,605,650,50),ISELCT 
@ 
C MTF OF THE DETECTOR IS THE FOURIER TRANSFORM OF ITS 
C SHAPE, A SINC FUNCTION. 
Cc 
602 WRITE (*,2) 
WRITE (*,2) 
WRITE (*,*) ‘ THIS TAKES TIME, SO BE PATIENT ’ 
$8(1)=1 
DO 603 I=2,257 
$8 (I) =SIN(3.14159*0.0003*18.65*I)/ 
x (3.14159*0.0003*18.65*I) 
TE ( Sone 0.2) SO(T) =o 
$8 (514-1) =S8(T) 
603 CONTINUE 
GO TO 609 
Cc 
C THIS SECTION BUILDS A GAUSSIAN APPROXIMATION TO THE MTF 
C CURVE FOR THE NPS-IRSTD X-DIRECTION OPTICS. 
Cc 
605 WRITE (*,2) 
WRITE (*,2) 
WRITE (*,*) ’ THIS TAKES TIME, SO BE PATIENT ’ 
$8(1)=1 
DO 606 I=2,257 
$8 (I) =EXP (-1.0* (I**2) /72000) 
Te (SO (T) . Li ee. 2)eS8 (1) =0 
$8 (514-1) =S8 (I) 
606 CONTINUE 
fe 


€ WITH THE MTFS CALCULATED, EACH ROW IS REMOVED FROM THE 
c ARRAY AND PROCESSED USING THE CONVOLUTION THEOREM BEFORE 
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c 
c 
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BEING PUT BACK INTO THE IMAGE. 


609 DO 610 J=0,449,5 
CALL ISGETP (NFRAME,J,0,512, LINE) 
DO 620 K=1,512 
RELINE (K) =LINE (K) -128.0 
620 CONTINUE 
N=512 
ID=1 
N1=9 
CALL PREPRO(RELINE,CG,N) 
CALL FE Tee, Ni) 


DO 625 L=1,512 
IF (S8(L) .LT. 0.3) THEN 
Ger) —0 
GO TO 625 
ENDIF 
GG(L)=cG (by / Ss) 
625 CONTINUE 
ID=0 
CALL FFT(CG,ID,N1) 
CALL POSPRO (CG, RELINE, AMAG,N) 
DO 630 M=1,512 
LINE (M) =INT (RELINE (M) +128) 
IF (LINE(M) .GT. 255) LINE(M)=255 
630 CONTINUE 
DO 635 L=0,4 
CALL ISPUTP (NFRAME, J+L, 0,512, LINE) 
635 CONTINUE 
610 CONTINUE 
CALL ISOTFR(NFRAME) 
GO TO 601 


THIS SECTION BUILDS A GAUSSIAN APPROXIMATION TO THE MTF 
CURVE SEOR, THE NES-IRStD YY DOrRBerLOneer Ties. 


650 WRITE (*,2) 
WRITE (*,2) 
WRITE (*,*) ’ THIS TAKES TIME, AND LOTS OF IT ’ 
WRITE (*,*) ‘ SO SIT BACK AND RELAX ’ 
WRITE (*,2) 
$8(1)=1 
DO 652 I=2,257 
S8 (I) =EXP (-1.0* (I**2) /3.27E+6) 
DE (SClGl ee OleZ) ee Siai( is) 0 
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$8 (514-1) =S8 (I) 


652 CONTINUE 


WITH THE MTF CALCULATED, EACH ROW IS REMOVED FROM THE 
ARRAY AND PROCESSED USING THE CONVOLUTION THEOREM BEFORE 
BEING PUT BACK INTO THE IMAGE. 


662 


664 


675 


677 


DO 660 J=1,512 


WRITE (*,/ (A18,2X,13,2X,A7)‘') ‘ WORKING ON COLUMN’, 
Weil,’ OW Bes" 
DO 662 K=0,449,5 
CALL ISGETP (NFRAME,K,J-1,1, LDET) 
RELINE ( (K/5) +1) =LDET-128.0 
CONTINUE 
DO 664 K=91,512 
RELINE (K) =0 
CONTINUE 
N=512 
i 
N1=9 
CALL PREPRO(RELINE,CG,N) 
CAD FET(CG, ID, N1) 
DO 675 L=1,512 
Te (SS (L)) .LT. 0.3) THEN 
CG (L) =0 
GO TO 675 
ENDIF 
CG (L) =CG(L) /S8 (L) 
CONTINUE 
ID=0 
CALL PERG mh Nil ) 
CALL POSPRO(CG, RELINE, AMAG, N) 
DO 677 K=0,449,5 
LDET=INT (128+RELINE ( (K/5) +1) ) 
Te (LDET .GT. 255) @aDET=255 
CALL ISPUTP (NFRAME,K, J-1,1,LDET) 
CALL ISPUTP (NFRAME, K+1,J-1,1,LDET) 
CALL ISPUTP (NFRAME, K+2,J-1,1,LDET) 
CALL ISPUTP (NFRAME, K+3,J-1,1, LDET) 
CALL ISPUTP (NFRAME, K+4,J-1,1,LDET) 
CONTINUE 


660 CONTINUE 


Sento 601 
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THIS SECTION SMOOTHS OUT THE 1 TO °SeROW EAPANCEChHeUSEr 
IN THE DMA PROCESS TO FILL UP THE DISPLAY. THE ROWS ARE 
WEIGHTED TO LIMIT THE OVERALL EFFECT. 


WRITE (*,*) ‘ENTER THE NUMBER OF THE FRAME TO SMOOTH ’ 
READ (*,’(I2)’) NFRAME 

WRITE (*,2) 

WRITE (*,2) 

WRITE (*,*) ’ THIS TAKES TIME, SO BE PATIENT ’ 


DISPLAY THE SELECTED FRAME 


CALL ISDISF (1) 
CALL ISOTFR (NFRAME) 


PULL 5 ROWS OUT, AND APPLY A WEIGHTED AVERAGE TO 
RECALCULATE THE VALUES FOR THE CENTER ROW. 


SELECT FRAME TO CORRECT. 


700 WRITE (*,*) ’ ENTER NUMBER OF FRAME TO CORRECT ’ 


705 


* 


READ (*,’(I2)’) NFRAME 

DOr Oil =27149 
CALL ISGETP (NFRAME,I~-2,0,512,IROW1) 
CALL ISGETP (NFRAME,I-1,0,512, IROW2) 
CALL ISGETP (NFRAME,I1,0,512, IROW3) 
CALL ISGETP (NFRAME,I+1,0,512, IROW4) 
CALL ISGETP (NFRAME, 1+2,0,512,IROWS) 
DO 705) u=0, 517 

TROW3 (J) =(IROW1 (J) +2* IROW2 (J) +3* IROW3 (J) + 
2* TROW4 (J) +IROWS (J))/9 

CONTINUE 
CALL ISPUTP (NFRAME,1,0,512, IROW3) 


701 CONTINUE 


GO TO 50 


THIS SECTION ALLOWS SAVING A FRAME BUFFER TO A FILE. 
IT AUTOMATICALLY ENCODES DATE, TIME, AND OTHER DATA. 
SEE APPENDIX C OF THE DT-IRIS MANUAL. 


800 WRITE (*,*) ‘ WHICH FRAME WOULD YOU LIKE TO SAVE ’ 


WRITE (*,*) ‘ ENTER VALUE 0 TO 15 ’ 

READ (*,’(1I2)’) NFRAME 

IF (NFRAME .LT. 0 .OR. NFRAME .GT. 15) GO TO 800 
WRITE (*,*) ’ TYPE THE NAME, INCLUDING EXTENSION ’ 
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WRITE (*,*) ° OF THE FILE TO WRITE THE FRAME ’ 
READ (*,’ (A30)’) FILES 
CALL ISSAVE (NFRAME,0,1,0,FILES) 


GO Toe 50 
ce 
€ THIS SECTION ALLOWS RESTORING DT-IRIS FILES TO A FRAME 
Cc TSUNI IME ISS 
Cc 
900 WRITE (*,*) ’ ENTER NAME OF FILE TO RETRIEVE ' 
READee, (220) 7 ) ETLES 
WhEEEe(*,*) “ RESTORE TO FRAME ? (VALUE © TO 15) 
READ (*,’ (I12)’) NFRAME 
IF (NFRAME .LT. 0 .OR. NFRAME .GT. 15) GO TO 900 
CALL ISREST (NFRAME,0,0,FILES) 
CALL ISOTFR (NFRAME) 
GG TO 50 
c 
C THIS SECTION ALLOWS ONE BUFFER TO BE COPIED TO A SECOND 
Cc ODT IT IER 
Cc 


LOCC e hikes, >) © ENTER® THE NUMBER @F THE SOURCE BUFFER ° 


READ (*,'’(I2)’) NFRAME 
IF (NFRAME .LT. 0 .OR. NFRAME .GT. 15) GO TO 1000 
WRITE (*,*) ’ ENTER THE NUMBER OF THE DESTINATION 
* BUFFER ’ 

READ (*,*) IDEST 

(Oeste. O .OR. IDEST .GT. 15) co Toe 1000 
CALL ISFCOP (NFRAME, IDEST) 

CALM ISDISP (1) 

CALL ISOTFR(IDEST) 


WRITE (*,*) ‘YOU ARE NOW LOOKING AT THE DESTINATION 


i BUFFER ’ 
GO" TOrs0 


ANOTHER. 


Ore t@n@) 


THIS SECTION ALLOWS ONE BUFFER TO BE SUBTRACTED FROM 


1TOO WRITE (x,*) ° ENTER THE NUMBER OF THE BUFFER TO BE 


* SUBTRACTED ' 
READ (*,’(I2)') NFRAME 
IF (NFRAME .LT. 0 .OR. NFRAME .GT. 15) GO TO 1000 


WRITE (*,*) ’ ENTER THE FRAME TO SUBTRACT IT FROM ’ 


READ (*,*) IDEST 

Th e(IDEST .LT. 0 .OR. IDEST .GT. 15) GO TO 1000 
CALL ISSUB(NFRAME, IDEST, IDEST) 

GAL 1SbiSr (1) 
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CALL ISOTFR (IDEST) 
WRITE (*,*) *’ YOU ARE NOW LOOKING BUFEER IDES. 


Cerro 50 
C 
Cc PROGRAM ENDS. ISEND MUST BE USED TO PREVENT POSSIBLE 
c CONTAMINATION OF SUBSEQUENT PROGRAMS. 
© 
1500 CALL ISEND 
STOP 
END 
C 
Cc SUBROUTINES IN THIS SECTION WERE WRITTEN BY PROF. E.A. 
Cc MILNE AND USED FOR THE LAB II CLASS (PH 2013) FOR 
c INSTRUCTION ON USE OF THE FAST FOURIER TRANSFORM. 
cS 
SUBROUTINE FFT(CG,ID,N1) 
C FOURIER TRANSFORM SUBROUTINE Latest update 3/14/88 
Cc by E. A. Milne 
c PARAMETERS IN THE CALL 
Se CG(*) = COMPLEX INPUT ARRAY AND OUTPUT ARRAY 
S ID = 1 FOR FORWARD TRANSFORM AND ID = 0 FOR INVERSE 
Cc TRANSFORM. 
Cc Nl IS A MEASURE OF THE SIZE OF THE ARRAY TO BE 
Sc TRANSFORMED. 
Ce SQ=SQRT (FLOAT (N2)) = NORMALIZING FACTOR 
Cc 
c LOCAL VARIABLES USED 
(Ss I, LL, M, kK, L, J, AND N = INTEGERS USED AS ARPAY 
c INDICES VOR, OR, DOC LUCE INDICES. One bo li 
c KL AND KU ARE INTEGERS FOR THE LOWER AND UPPER LIMITS OF 
c A DO LOOP. 
Cc N2 = 2**N1 = THE TOTAL NUMBER OF TERMS IN THE ARRAY. 
c TH = 2*PI/N2 = PHASE ANGLE INCREMENT. 
iS 
c TYPE DECLARATION 
iS 
IMPLICIT COMPLEX (C) 
(MPEG UE INTEGER*2 (i-N)) 
DIMENSION CG(1) 
eC 
Sc N2 AND TH ARE CALCULATED 
¢ 


N2=2**N1 
TH=6 .283185/N2 
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THIS SECTION IS THE SUMMING SECTION - Nl SUMS 


2 DO Ss M=1,N1 
h=2** (1 -—M) 
L=2** (M-1) 
DO 4 LL=1,L 
ne he 
Clr SrETINY (KC,N, Nl-1) 
Cb=CExE (CMP ELA (Ol, TA* FLOAT (N))) 
IF (ID.EQ.1) CD=CONJG(CD) 
hae * KC*Kt 
KU=KL+K-1 
DO 3 I=KL, KU 
CA=CG (I) +CG (I+K) *CD 
Cengtik) —CG (1) -CG (1K) *CbD 
CG(I)=CA 
3 CONTINUE 
* CONTINUE 
5 CONTINUE 


THE REORDERING AND NORMALIZATION IS DONE IN THIS SECTION 


IF (ID.EQ.1) THEN 
SQ=1.0 

BEGE 
SQ=F LOAT (N2) 

ENDIF 

DO 8 I=1,N2 
J=I-1 
CALL BITINY (dc, N,N1) 
ECON) Onc 


6 N=N+1 

CA=CG (I) /CMPLX (SQ, 0.0) 

CG (I) =CG (N) 

CG (N) =CA 
ai CG (I) =CG (I) /CMPLX(SQ,0.0) 
8 CONTINUE 


NOW RETURN 


RETURN 
END 
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SUBROUTINE BITINV(K,L,M) 
FINDS THE BIT INVERSION OF A (M) BIT NUMBER (K) TO (L) 
THIS IS USED IN THE FOURIER TRANSFORM SUBROUTINE 

IMPLICIT INTEGER*2 (I-N) 

J=K 

L=0 

DO 1 I=1,M 

L=L*2+J-J3/2*2 

1 ded/ 2 

RETURN 

END 


SUBROUTINE PREPRO(RE,CG,N) 


THIS SUBROUTINE CONVERTS A REAL ARRAY RE TO A COMPLEX 
ARRAY CG TO BE USED BY THES PPT ROULINE. 


IMPLICIT? SN PEGERs 2 b-N) 
IMP LIC UPR RCOMP EEX MC) 
DIMENSION RE(N) ,CG(N) 
bes i= 1, N 

CG (I) =CMPLX (RE (I), 0.0) 

1 CONTINUE 

RETURN 
END 


SUBROUTINE POSPRO (CG, RE, AMAG, N) 


THIS SUBROUTINE TAKES A COMPLEX ARRAY AND FINDS THE REALS 
RE, THE IMAGINARIES AIM, THE MAGNITUDES AMAG, THE PHASES 
PH, AND THE POWER SPECTRUM PO@S@—1HAT THEY -CAN OBE 
SEPARATEEY PEOTIEDR: 


IMPLICIT INTEGER*2 (I-N) 
IMPLICIT COMPLEX (C) 
DIMENSION CG(N),RE(N) , AMAG(N) 
DO lt 1=1, N 

RE (I)= REAL (CG(I) ) 

AMAG (I) =CABS (CG (I) ) 

1 CONTINUE 

RETURN 
END 
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SUBROUTINE HAMM(RS,RD,N1) 


THES SUBROUTINE PARTIALLY REMOVES THE DISCONTINUITY 
BETWEEN THE FIRST AND LAST POINTS OF THE ARRAY RS AND 
STORES THE RESULT IN THE ARRAY RD. AS A NOTE TO FUTURE 
USERS OF THIS SYSTEM WHEN THIS ROUTINE WAS INCLUDED IN 
Bie PRECESSING, IT PRODUCED SEVERE INTENSITY REDUCTIONS 
IN THE SIDES OF THE IMAGE. 


THE NUMBER OF POINTS IN THE ARRAY = N = 2**N1 


IMPLICIT INTEGER*2 (I-N) 
DIMENSION RS(1),RD(1) 
N=2**N1 
DOO I = 1,N 
RD (I) =RS (I) * (0.54-.46*COS (6.283185*FLOAT (I-1) / 
* FLOAT (N-1) ) ) 
10 CONTINUE 
RETURN 
END 


SUBROUTINE HANN (RS, RD,N1) 


THIS SUBROUTINE COMPLETELY REMOVES THE DISCONTINUITY 
BETWEEN THE FIRST AND LAST POINTS OF THE ARRAY RS AND 
STORES THE RESULT IN RD. 


THE NUMBER OF POINTS IN THE ARRAY = N = 2**N1 


IMPLICIT INTEGER*2 (I-N) 
DIMENSION RS(1),RD(1) 
N=2**N1 
DOmOnt = 1,N 
RD (I) =RS (I) *0.5* (1.-COS (6.283185*FLOAT (I-1)/ 
* FLOAT (N-1)) ) 
10 CONTINUE 
RETURN 
END 
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APPENDIX D 


I. NPS-IRSTD DETECTOR CALIBRATION DATA 
A. RELATIVE RESPONSIVITY RESULTS 


Side 1 Amplitude (volts) Side 2 Amplitude (volts) 
101 822 3 BOIL 0.020 
102 I 5 SSO 202 3.433 
103 IL, Stake, 203 35 382 
104 1.904 204 35 18S 
105 1.900 205 35 Oz 
106 iL Pay 206 3.400 
iLO) 7 1.681 207 3, 26d 

* 108 0.071 208 Shp Sia al 
109 io @Hil 209 3.302 
110 ie oZe 210 3.163 
ial SL oe 211 So lSzZ 
iL LZ iL 6 OSL 212 3.320 
llas 1 SS Als} 35645 
114 sce) 214 3, es 

<r lelesy 0.400 215 3, 180 
116 A tesyk ZanG 3.166 
ill ae? IL BSN ple) 3,262 
SL AL} i, VEZ 218 3, 30S 
119 er Oe 219 3239 
LAC 1.691 220 3 Lis: 
Lei IL. ASO BAI 2.802 
122 1.763 222 3.028 
123 1.612 BAS 8), 202 
124 Lo FOAL 224 Sree) 
EAS 1.643 225 3.140 
126 1,992 x 226 OROZG 
Ley 1.642 * 227 O.022 
128 1.661 * 228 0.600 
129 iL, SAL * 229 0.600 

we ALSO 0.050 * 230 0.030 
iL Shai i 4 Gib aL if 2 Sh IL 3.050 
IL Sz 1.692 Oe 24S 
SS iL 5 FOZ Bs 3.050 
iL shal i, HAZ 234 3.000 


* Denotes that the detector channel’s response is 
irregular. 
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Side 1 Amplitude (eles) SHiclionmne Amplitude (volts) 


135 600 230 Sse 
13'S os 0 236 S120 
MEG oo 0 237 S122 
138 ice l 238 3.030 
139 roo 0 239 3150 
140 1.740 240 Sepsis Il 
141 1.604 241 Seco 
142 1 730 242 3.410 
143 16 243 3.450 
144 te Biles 244 3.452 
145 eos 245 3.404 
146 iene eae 246 3.303 
147 1.840 247 3.233 
148 Ths alles 248 See or 
149 IL te sal 249 Seoes 
Asie) (DEEN) S10) PD) dee 3 
aya i oGa Bis ll Seow 
a2 oo 9 292 322068 
L383 v.71 o4 203 Sells 
154 83.5 254 Ss. 201 
oo ieee 2090 Cee a 
tS‘ dr 256 3.144 
ey) ino Zod See) ) 
Jey! Ay TO 200 3.280 
ES, 17 O.1 209 Sor 
160 ee 260 Seo! 
Low 16 90. 261 3.058 
62 I ais} 262 Bos b AL 
nes ied 263 S102 
164 Ie Pee 264 S055 
Jey) 1.742 265 3.141 
166 vo Le 266 5) Oe IL 
LG? Nl esits) 267 oral Od 
168 Te733 268 S42 
Wee) 820 2Go e267 
Ig) i 20 270 ery illo) 
171 7 70 eee al GO. 015 
72 er oe. ae 3.340 
dhs) 1.824 213 Bec 9 
174 ess 274 So 2S) 
das UEOse 215 3.160 
dL IS 0.030 276 3.042 
dea a 829 Ey 3.164 
i738 W5030 27s 2713 
iS To 02 7719 22080 
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Side_1 Amplitude (volts) Side 2 Amplitude (volts) 


180 T2270 ~ 230 2.480 
181 1.676 261 3.142 
182 1.854 282 2.900 
183 1.704 265 2.882 
184 ereswanl 284 S.052 
185 Ale ie 285 Sez 
186 E.7St 286 3. 330 
187 £2662 Patel) So 16a 
188 T2662 288 22952 
189 1.640 239 2.840 
190 1,681 290 22019 


* Denotes that the detector channel’s response is 
irregular 


B. LIST OF DEFECTIVE DETECTOR CHANNELS 


Lag Array Lead Array 
O01 08 
26 15 
BY 30 
28 50 
29 25 
30 76 
50 78 
im 

x 78 
719 
80 


x Denotes that detector is operational but signal strength 
is significantly diminished. 
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